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INTRODUCTION 

This  manual   contains   input   instructions,    recommendations,   and 
illustrations   for  using  a  new  culvert  analysis /design  methodology, 
called  CANDE    (Culvert  ANalysis   and  DEsign) .      CANDE  is   a  user-oriented, 
FORTRAN  computer  program  employing  state-of-the-art   analytical 
methods   for  designing  and  analyzing  the  structural  integrity  of  buried 
pipe  culverts. 

For  convenience  and  clarity,   this  manual  Is   divided  into  three 
parts.     Part  I  is  a  brief  summary  of  the  scope  of  the  CANDE  program. 
Part   II  provides   formated  input  instructions   for  data  cards   along  with 
extensive   commentaries   and  recommendations.      Part   III  contains  example 
problems  covering  a  wide  range  of  potential  applications  and  illustrating 
input /output   features   of  CANDE.      For  detailed  developments   the   reader 
is  referred  to  the   companion  docun\    ':s: 

•  CANDE:    Engineering  Manua.         A  Modern  Approach  for  the   Structural 
Design  and  Analysis   of  E  d  Cylinders,   by  M.    G.    Katona  et   al. 

Civil  Engineering  Labors    ,-ry,  Port  Hueneme,   CA, 1976 

(final  technical  report   lo  Federal  Highway  Administration) 

•  CANDE:    System  Manual,   by  M.    G.    Katona  and  J.   M.    Smith.    Port 
Hueneme,   CA, 1976    (Report  to  Federal  Highway  Administration) 


PART  I 
CANDE  OVERVIEW 

CANDE  offers  a  choice  of  three  solution  levels  together  with  a 
choice  of  five  pipe  types  for  designing  and/or  analyzing  buried  pipe 
culverts.  The  solution  levels  range  from  an  elasticity  solution  to 
a  general  finite  element  solution,  which  enables  the  engineer  to  select 
a  degree  of  rigor  and  cost  commensurate  with  the  project  at  hand.  The 
pipe  types  include  corrugated  steel,  corrugated  aluminum,  reinforced 
concrete,  plastic  pipe,  and  basic  (nonstandard). 

Figure  1-1  presents  a  schematic  overview  of  CANDE  that  shows  the 
three  principal  areas:  the  main  control  area  at  top,  the  pipe  library 
at  bottom  left,  and  the  solution  library  at  bottom  right.  The  main 
control  area  reads  the  problem  control  input  and  then  acts  as  a  switch- 
board by  shifting  information  back  and  forth  between  the  pipe  library 
and  the  solution  library. 

To  start  each  problem,  three  fundamental  selections  are  required 
by  the  user:   (1)  execution  mode,  (2)  solution  level,  and  (3)  pipe 
type. 

Execution  mode  is  the  decision  between  design  or  analysis.  Analysis 
means  a  particular  pipe  soil  system  is  completely  defined  and  then 
solved  by  the  chosen  solution  level.  Output  consists  of  the  structural 
responses  (displacements,  stresses,  strains)  as  well  as  an  evaluation 
of  the  pipe  performance  in  terms  of  factors -of- safety  against  potential 
modes  of  failure. 


CANDE...CuIvert  Analysis/ Design 
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Figure  1-1.  Schematic  diagram  for  CANDE. 


The  alternative  execution  mode,  design,  requires  the  same  input 
definition  except  the  pipe  wall  geometrical  section  properties  will  be 
unknown.   Instead,  desired  safety  factors  are  input,  and  CANDE  achieves 
a  design  by  a  direct  search  approach.  That  is,  a  series  of  analyses 
are  performed  such  that  an  initial  trial  section  is  successively  modified 
until  the  desired  safety  factors  are  achieved.  Design  output  includes 
required  wall  properties,  actual  factors  of  safety,  and  structural 
responses.  Naturally,  the  required  wall  properties  depend  on  the  pipe 
type.  For  example,  the  properties  for  metal  pipes  are  given  in  corrugation 
and  gage  sizes,  while  the  reinforced  concrete  pipe  properties  are  given 
in  wall  thickness  and  steel  area. 

The  three  solution  levels  in  the  solution  library  correspond  to 
sucessively  increased  levels  of  analytical  power  and  sophistication. 
The  successive  increase  in  analytical  power  is  accompanied  by  an  increase 
in  solution  costs.  All  solution  levels  assume  plane  strain  geometry 
and  are  cast  in  incremental  form  to  simulate  incremental  construction. 

Level  1  is  an  enhanced  elasticity  solution  applicable  for  round 
pipes  deeply  buried  in  homogeneous  soil.  Levels  2  and  3  share  a  common 
finite  element  methodology  and  only  differ  on  the  manner  of  input  - 
automatic  or  manual.  Level  2  provides  a  "canned*'  finite  element 
mesh  applicable  for  most  symmetric  culvert  installations,  including 
trench  and  embankment  conditions.  Level  3  permits  consideration  of 
any  arbitrary  culvert  configurations;  however,  the  finite  element  mesh 
model  must  be  defined  by  the  user. 


Table  1-1  summarizes  the  analytical  features  and  limitations  of 
each  solution  level  along  with  an  indication  of  input  requirements  and 
computer  running  time . 

The  CANDE  pipe  library  currently  contains  subroutines  for 
corrugated  steel,  corrugated  aluminum,  reinforced  concrete,  and  smooth 
wall  plastic  pipe.   In  addition  a  subroutine  called  * 'BASIC* *  allows 
for  description  of  nonstandard  pipes;  however,  this  routine  is  only 
applicable  to  analysis  problems.   The  pipe  subroutines  service  all 
three  solution  levels  and  are  the  key  control  areas  of  CANDE.  They  perform 
four  operations:   (1)  read  pipe  input  data,  (2)  modify  pipe  properties 
for  nonlinear  responses,  (3)  evaluate  pipe  safety  factors,  and  if  in 
design  mode,  (A)  update  section  properties  of  pipe  wall.  Table  1-2 
summarizes  the  basic  characteristics  of  each  pipe  subroutine. 

The  foregoing  view  of  the  CANDE  program  illustrates  the  "mix- 
and-match"  features  of  CANDE.  That  is,  any  pipe  type  from  the  pipe 
library  can  be  matched  with  any  solution  level  from  the  solution  library, 
and  the  pair  can  be  run  in  either  a  design  or  analysis  mode.  Little 
or  nothing  has  been  said  about  the  various  technical  options,  such 
as,  incremental  construction,  nonlinear  soil  laws,  friction  interface 
models,  nonlinear  pipe  models,  and  other  modeling  features.  These  capa- 
bilities will  be  discussed  in  the  appropriate  sections  of  the  user 
instructions  (Part  II) ,  and  they  are  discussed  in  depth  in  the  engineering 
manual . 
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Table  1-2.   Characteristics  of  Pipe 


Type 
of 

Material 
Stress-Strain 

Design/Evaluation 

Design 

Pipe 

Models 

Criteria 

Output 

Corrugated 

Linear;  yield- 

Excessive  deflection; 

List  of  required 

steel 

hinge  theory; 

thrust  yielding; 

gage  thickness 

bilinear  stress- 

elastic  buckling*2 

for  each  standard 

strain 

steel  corrugation 
size 

Corrugated 

Linear;  yield- 

Excessive  deflection; 

List  of  required 

aluminum 

hinge  theory; 

thrust  yielding; 

gage  thickness 

bilinear  stress- 

elastic  buckling*3; 

for  each  standard 

strain 

outer  fiber  strain 

aluminum  corruga- 

rupture 

tion  size 

Reinforced 

Concrete 

Concrete  crushing; 

Required  steel 

concrete 

cracking  and 

steel  yielding; 

area  for  double 

nonlinear  com- 

concrete shear 

circular  or 

pression;  yield- 

stress; crackwidth, 

elliptical  rein- 

ing of  steel 

0.01  in. 

forcement  cages; 

reinforcement 

required  wall 
thickness 

Smooth 

Linear 

Excessive  deflection; 

Wall  thickness 

wall 

outer  fiber  stress; 

plastic 

elastic  buckling 

Basic 

Linear 

No  design, 

No  design, 

analysis  only 

analysis  only 

Buckling  prediction  is  an  approximation  external  to 
the  solution  levels. 


PART  II 
INPUT  INSTRUCTIONS  AND  COMMENTARIES 

For  each  problem,  input  data  must  be  prepared  in  accordance  with 
instructions   from  three  separate  sections.     The  purpose  of  each  section 
is: 

Section  A  -  Master  control. 

Section  B  -  Pipe-type  definition.     This  section  is  subdivided  into: 
corrugated  steel,   corrugated  aluminum,   reinforced 
concrete,  plastic,  and  basic. 
Section  C  -  Solution  level  system  definition.     This  section  is 

subdivided  into:   Level   1    (modified  elasticity  solution), 
Level  2    (finite  element  solution,   automated),   and 
Level  3   (finite  element  solution,   standard). 

These  sections  are  presented  in  order  with  the  subsections  identified 
by  name.      Data  card  input  instructions  are  grouped  for  each  subsection, 
followed  by  a  narrative  commentary.   The  input  data  are  keyed  to  the 
commentary  by  numbers  in  the  right-hand  column.     Each  problem  must  begin 
at   Secti*n  A  and  follow  the  instructions  through  Section  C. 


SECTION  A  -  MASTER  CONTROL  CARD 

Input 

Card  1A.   Master  control  card  (one  card  per  problem) 


Columns 

Variable 

(format) 

(units) 

01-06 

XMODE 

(A4,2X) 

(word) 

Entry  Description 

Word  defining  program  mode, 

=  ANALYS,  denotes  analysis  problem 
=  DESIGN,  denotes  design  problem 
=  STOP,  program  terminates,  last 
card  in  deck 


Note 
(1) 


08-08 
(ID 


LEVEL 


Defines  solution  level  to  be  used, 
=  1,  denotes  elasticity  solution 
=  2,  denotes  finite  element 

solution  with  automated  mesh 
=  3,  denotes  finite  element 

solution  with  user-defined  mesh. 


(2) 


10-15 
(A4,2X) 


PTYPE 
(word) 


17-76 

HED 

(15A4) 

(words) 

77-78 

NPMAT 

(12) 

79-80 

NPPT 

(12) 

Defines  pipe  type  to  be  used,  (3) 

=  STEEL,  denotes  corrugated  steel 
=  ALUMIN,  denotes  corrugated 

aluminum 
=  CONCRE,  denotes  reinforced 

concrete 
=  PLASTI,  denotes  smooth  plastic 
=  BASIC,  denotes  arbitrary  pipe, 

for  analysis  only 

User  defined  heading  of  problem  to  be 
printed  with  output 

Number  of  pipe  elements;  only  required 
when  LEVEL  =  3 

Number  of  pipe  nodes;  only  required  when 
LEVEL  =  3 


***  GO  TO  SECTION  B  *** 


Commentary 

(1)  The  decision  of  design  or  analysis  is  controlled  by  specifying 
the  variable  XMODE.     Analysis  implies   all  system  and  pipe  properties 

are  known,   and  the  objective  is  to  evaluate  the  pipe  performance.     Design 
implies  the  pipe  wall  section  properties   are  unknown,   and  that  they  will 
be  the  object  of  the  design  search.     Note,   the  program  will  continue 
to  execute  problems  back-to-back  until  XMODE  «■  STOP  is  encountered. 

(2)  Selection  of  a  solution  level  is   a  question  of  what  degree 
of  detail  is  desired  and  how  much  cost  and  time  is  it  worth.   Level  1 
is  quick  and  comparatively  simple,  whereas  Level  3  is  detailed  and 
sophisticated.      Table   1-1   of  Part  I  summarizes  the  scope  of  each  solution 
level,   and  Section  C  of  Part  II  provides  additional  detail.     As  a  general 
rule  of  thumb,  one  should  use  the  solution  level  commensurate  with 

the  confidence  in  assumed  values   for  loading,   soil  properties,   geometry 
and  other  culvert  aspects. 

For  many  common  design  problems,  Level   1   is   commensurate  with 
knowledge  of  the  pipe-soil  system  and  is   a  good  design  tool.   However, 
for  shallow  burial  depths,  i.e.,   less   than  one  radii  of  cover,  Level 
1   should  not  be  used. 

Level  2,  which  is  considered  the  workhorse  of  CANDE,  is  applicable 
to  most  design  and  analysis  problems.      Input  is  relatively  simple,  yet 
a  great  deal  of  modeling  detail  is  possible. 
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A. special  * 'extended' *  feature  of  Level  2  allows  selective  modification 
of  the  standard  Level  2  mesh. 

Level  3  requires  familiarity  with  the  finite  element  method  for 
proper  modeling  of  the  pipe-soil  system.  This  level  is  required  for 
nonstandard  problems,  such  as  asymmetric  loading  and  geometry.  Level 
3  unleashes  the  full  power  of  the  finite  element  method;  however, 
a  considerable  amount  of  time  is  required  to  define  input  and  "debug** 
the  mesh. 

(3)  Table  1-2  of  Part  I  summarizes  characteristics  of  each  pipe  type. 
Only  one  pipe  type  can  be  selected  for  each  problem.  To  obtain  a  com- 
parative study  of  one  pipe  type  versus  another,  simply  resubmit  the 
data  deck  with  the  new  PTYPE  and  corresponding  Section  B  data  cards. 

If  no  pipe  elements  are  desired  (only  possible  in  Level  3),  PTYPE 
must  be  defined  as  BASIC,  and  NPMAT  and  NPPT  must  be  set  to  zero. 
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SECTION  B  -  PIPE  TYPES 

Go  to  the  subsection  that  corresponds  to  the  chosen  pipe  type: 

©  Corrugated  steel 

©  Corrugated  aluminum 

©  Reinforced  concrete 

©  Plastic 

©  Basic 

Commentaries  are  included  after  the  input  instructions  for  each  pipe 
type. 
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Corrugated  Steel  Input  (Cards  IB,  2B) 

Card  IB.   Pipe  size  and  material  properties: 


Columns 

Variable 

(format) 

(units) 

01-05 

NONLIN 

(15) 

06-15 

PDIA 

(F10.0) 

(in.) 

16-25 

PE 

(F10.0) 

(psi) 

26-35 

PNU 

(F10.0) 

36-45 

PYIELD 

(F10.0) 

(psi) 

46-55 

PDEN 

(F10.0) 

(pci) 

56-65 

PE2 

(F10.0) 

(Rsi) 

Entry  Description  Note 

Denotes  type  of  material  behavior       (1) 
=  0,  linear  stress  strain-law 
=  1,  yield-hinge  theory 
=  2,  bilinear  stress-strain  law 

Pipe  diameter,  average  (2) 


Elastic  Young's  modulus  of  pipe 
material, 

Default  =  30  x  10  psi 

Poisson's  ratio  of  pipe  material, 
Default  =0.3 

Yield  stress  of  pipe  material, 
Default  =  33,000  psi 

Density  of  material  for  body  weight; 
only  applies  to  Levels  2  and  3, 
Default  =0.0 

Modulus  of  upper  portion  of  (3) 

bilinear  stress-strain  curve;  used 
only  with  NONLIN  =  2 
Default  =0.0 
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C^rd  2B.  For  analysis  only,  wall  properties: 

Columns  Variable 

(format)  (units)      Entry  Description                     Note 

01-10  PA  9        Area  of  pipe  wall  section              (4) 

(F10.0)  (in.  /in.)   per  unit  length 

11-20  PI   ,        Moment  of  inertia  of  pipe 

(F10.0)  (in.  /in.)   wall  section  per  unit  length 

21-30  PS          Section  modulus  of  pipe  wall 

(F10.0)  (in.  /in.)   per  unit  length 

Card  2B.  For  design  only,  desired  safety  factors: 


Entry  Description  Note 

Desired  safety  factor  against  (5) 

thrust  yielding, 
Default  =3.0 

Desired  safety  factor  against 
20%  deflection  (20%  of  diameter), 
Default  =4.0 

Desired  safety  factor  against 
elastic  buckling, 
Default  =2.0 


Columns 

Variable 

(format) 

(units) 

01-10 

PFS(l) 

(F10.0) 

11-20 

PFS(2) 

(F10.0) 

21-30 

PFS(3) 

(F10.0) 

***  GO  TO  SECTION  C  *** 


14 


Corrugated  Steel  Commentary 

(1)  The  parameter  NONLIN  indicates  the  degree  of  material  nonlinear! ty 
to  be  considered.  Most  metal  culverts  experience  outer  fiber  yielding 

due  to  bending  under  normal  loading  conditions.  Consequently,  the  linear 
approximation  (NONLIN  ■  0)  is  questionable  except  for  parameter  studies. 
As  a  general  rule,  the  bilinear  assumption  (NONLIN  ■  2)  should  be 
used.  However,  for  Level  1  solutions,  the  yield-hinge  theory  (NONLIN 
-  1)  is  preferable  if  more  than  50%  of  the  wall  section  yields.  This 
is  because  Level  1  uses  a  smeared  average  of  the  bilinear  modulus  that 
does  not  properly  represent  localized  yielding.  The  yield-hinge  theory 
was  developed  expressly  to  treat  local  yielding  (hinging)  for  Level 
1. 

(2)  For  round  pipe,  PDIA  is  the  pipe  diameter  measured  from  wall 
center  to  wall  center.  For  elliptical  pipe  (Level  2),  PDIA  is  the 
average  of  the  major  and  minor  diameters.   For  other  culvert  shapes 
(Level  3),  PDIA  is  a  nominal  average  diameter;  further  shape  definition 
is  required  in  Section  C. 

(3)  The  bilinear  stress-strain  model  and  associated  parameters 
are  shown  in  Figure  II- 1.   The  default  values  provided  for  the  material 
parameters  are  recommended  for  standard  corrugated  steel  culverts. 
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(4)  Only  include  this  card  if  XMODE  »  ANALYS.  Note,  all  units 
are  inches. 

(5)  Only  include  this  card  if  XMODE  -  DESIGN.   Normal  ranges 
for  specified  safety  factors  are: 

Thrust  yielding,  2.0  -  3.0 
Excessive  deflection,  3.5  -  4.0 
Elastic  buckling,  2.0  -  3.0 


Strain 


Figure  II-l.    Bilinear  stress-strain  parameters  for  steel. 
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Corrugated  Aluminum  Input  (Cards  IB,  2B) 
Card  IB.   Pipe  size  and  material  properties: 


Columns 

Variable 

(format) 

(units) 

01-05 

NONLIN 

(15) 

06-15 

PDIA 

(F10.0) 

(in.) 

16-25 

PE 

(F10.0) 

(psi) 

26-35 

PNU 

(F10.0) 

36-45 

PYIELD 

(F10.0) 

(psi) 

46-55 

PELONG 

(F10.0) 

(in. /in. ) 

56-65 

PDEN 

(F10.0) 

(pci) 

66-76 

PE2 

(F10.0) 

(psi) 

Entry  Description  Note 

Denotes  type  of  material  behavior,       (1) 
=  0,  linear  stress-strain  law 
=  1,  yield-hinge  theory 
=  2,  bilinear  stress-strain  law 

Pipe  diameter,  average  (2) 


Elastic  Young's  modulus  of  pipe 
materials, 

Default  =  10.2  x  10  psi 

Poisson's  ratio  of  pipe  material, 
Default  =  0.33 

Yield  stress  of  pipe  material, 
Default  =  24,000 

Strain  at  material  rupture,  (3) 

Default  =  0.05  in. /in. 

Density  of  material  for  body  weight; 
only  applies  to  Levels  2  and  3, 
Default  =0.0 

Modulus  of  upper  portion  of 
bilinear  stress-strain  curve;  used 
only  with  NONLIN  =  2, 
Default  =  600,000  psi 
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Card  2B.   For  analysis  only,  wall  properties: 


Columns 
(format) 

Olr-10 
(F10.0) 

11-20 
(F10.0) 

21-30 
(F10.0) 


Variable 
(units) 

PA  2 
(in.  /in.) 

PI  4 
(in.  /in.) 

PS   „ 


Entry  Description 

Area  of  pipe  wall  section 
per  unit  length 

Moment  of  inertia  of  pipe 
wall  section  per  unit  length 

Section  modulus  of  pipe  wall 
per  unit  length 


(in.  /in.) 
Card  2B.   For  design  only,  desired  safety  factors: 


Columns 
(format) 

01-10 
(F10.0) 


11-20 
(F10.0) 


21-30 
(F10.0) 


31-40 
(F10.0) 


Variables 
(units) 

PFS(l) 


PFS(2) 


PFS(3) 


PFS(4) 


Entry  Description 

Desired  safety  factor  against 
thrust  yielding, 
Default  =3.0 

Desired  safety  factor  against 
20%  deflection  (20%  of  diameter), 
Default  =4.0 

Desired  safety  factor  against 
elastic  buckling, 
Default  =2.0 

Desired  safety  factor  against 
strain  rupture, 
Default  =2.0 


Note 
(4) 


Note 
(5) 


***  GO  TO  SECTION  C  *** 
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Corrugated  Aluminum  Commentary 

(1)  Recommendations  for  specifying  NONLIN  are  the  same  as  for 
corrugated  steel  (comment  1).  In  general,  use  NONLIN  »  2. 

(2)  For  round  pipe,  PDIA  is  the  pipe  diameter  measured  from  wall 
center  to  wall  center.  For  elliptical  pipe  (Level  2) ,  PDIA  is  the 
average  of  the  major  and  minor  diameters.  For  other  culvert  shapes 
(Level  3) ,  PDIA  is  a  nominal  average  diameter;  further  shape  definition 
is  required  in  Section  C. 

(3)  The  concern  for  strain  rupture  is  the  primary  difference 
between  the  treatment  of  corrugated  aluminum  and  corrugated  steel. 

The  ductile  range  of  aluminum  terminates  in  strain  rupture.  Accordingly, 
excess  bending  strain  is  considered  unsafe  for  aluminum.   The  bilinear 
stress-strain  parameters  are  illustrated  in  Figure  II-2  along  with  the 
rupture  strain  parameter.  The  default  values  of  these  parameters  are 
recommended  for  standard  corrugated  aluminum  pipe* 

(4)  Only  include  this  card  if  XMODE  -  ANALYS.  Note,  all  units 
are  inches. 

(5)  Only  include  this  card  if  XMODE  -  DESIGN.   Normal  ranges  for 
specified  safety  factors  are: 
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Thrust  yielding,   2.0  •    3.0 
Excessive  deflection,   3.5   -  4.0 
Elastic  buckling,  2.0   -  3.0 
Outer  fiber  rupture,  2.0  -  3.0 


i 

1 

u 

1               ■ 

rupture 

--— 1 

j 

PYIELD 

r 

PE2_L- — — " 

1 

vt 

O 

</5 

/l 

i 

PE 

1 
1 

1 
1 

£v_ 

Strain 


PELONG 


Figure  \\-2.    Bilinear  stress-strain  parameters  for  aluminum. 
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Reinforced  Concrete  Inpu' 
Card  IB.   Pipe  size,  shar 


:rds  IB,  2B,  3B) 
and  nonlinear  character: 


Columns 

Variable 

(format) 

(units) 

01-10 

PDIA 

(F10.0) 

(in.) 

11-20 

PT 

(F10.0) 

(in.) 

22-25 

RSHAPE 

(A4) 

(word) 

26-30 
(15) 


31-40 
(F10.0) 


41-50 
(F10.0) 


51-60 
(F10.0) 


NONLIN 


STNMAT(l) 
in. /in. 


STNMAT(2) 
in . / in . 


STNMAT(3) 
in . / in . 


Entry  Description  Note 

Inside  pipe  diameter,  average  (1) 


Concrete  wall  thickness,  for  design 
this  is  left  blank  unless  a  fixed-wall 
design  is  desired 

Name  to  select  shape  of  steel  cage(s),   (2) 
=  CIRC,  circular  cages,  inner  and 

Outer 
=  ELLI,  elliptical  reinforcement, 

one  cage 
=  ARBI,  arbitrarily  shaped  cages, 
this  applies  to  analysis  only, 
Default  =  CIRC 


Degree  of  nonlinearity, 

=  1,  concrete  cracking  only 
=  2,  also  include  nonlinear 

compression  of  concrete 
=  3,  also  include  steel  yielding, 
Default  =  1 

Concrete  strain  at  which  tensile 
cracking  occurs  (positive) , 
Default  =0.0  in. /in. 

Concrete  strain  at  elastic  limit 
in  compression  (positive), 

Default  =  1/2  PFPC/PCE  (see  next 

card) 

Concrete  strain  at  initial 
compressive  strength,  f^  (positive), 
Default  =  0.002  in. /in. 


(3) 


(4) 
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Card  2B.   Concrete  and  steel  properties: 


Columns 
(format) 

01-10 
(F10.0) 


11-20 
(F10.0) 

21-30 
(F10.0) 

31-40 
(F10.0) 


Variable 
(units) 

PFPC 
(psi) 


PCE 
(psi) 

PNU 


PDEN 
(pcf) 


41-50 

PFSY 

(F10.0) 

(psi) 

51-60 

PSE 

(F10.0) 

(psi) 

61-70 

PSNU 

(F10.0) 

Enfry  Description 

Compressive  strength  of 
concrete,  f^., 

Default  =  4,000  psi 


Young's  modulus  of  concrete  in 
elastic  range,  ,      . 

Default  =  33 (density)  '(f1) 

Poisson  ratio  of  concrete, 
Default  =  0.17 

Unit  weight  of  concrete  (density), 
Default  =0.0  for  body  weight; 
however , 

Default  =  150  pcf  for  modulus 
calculation 

Yield  stress  of  reinforcing  steel, 
Default  =  40,000  psi 

Young's  modulus  of  steel 
Default  =  29  x  106  psi 

Poisson' s  ratio  of  steel, 
Default  =0.3 


Note 
(5) 
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Card  3B.   For  analysis  only,  steel  area  and  location: 


Entry  Description  Note 

Area  of  inner  cage  reinforcement        (6) 
steel;  If  RSHAPE  =  ELLI,  this  is 
area  of  the  elliptical  case 

Area  of  outer  cage  reinforcement 
steel;  not  used  if  RSHAPE  =  ELLI 

Thickness  of  concrete  cover  to 
center  of  inner  cage,  or  minimum 
cover  of  elliptical  cage, 
Default  =  1.25  in. 

31-40       TBO  Thickness  of  concrete  cover  to 

(F10.0)      (in.)        center  of  outer  case,  not  used  if 

RSHAPE  =  ELLI, 

Default  =  1.25  in. 


Columns 
(format) 

Variable 
(units) 

01-10 
(F10.0) 

ASI  2 
(in.  /in.) 

11-20 
(F10.0) 

ASO  2 
(in.  /in.) 

21-30 
(F10.0) 

TBI 
(in.) 
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48-50 
(A3) 


51-60 
(F10.0) 


61-70 
(F10.0) 


For  design  only,  design  specifications 


Entry  Description 


Card  3B. 

For  design 

Columns 

Variable 

(format) 

(units) 

01-10 

PFS(l) 

(F10.0) 

11-20 

PFS(2) 

(F10.0) 

21-30 

PFS(3) 

(F10.0) 

31-40 

ALCW 

(F10.0) 

(in.) 

GOAL 
(word) 


TBI 
(in.) 


SRATIO 


Desired  safety  factor  against 
steel  yielding, 
Default  =1.6 

Desired  safety  factor  against 
concrete  crushing, 
Default  =2.0 

Desired  safety  factor  against 
concrete  shear  failure, 
Default  =2.0 

Allowable  maximum  crack  width 
for  the  design  load;  normally 
taken  as  0.01  in, 

Default  =  Crack  width  is  not 
considered  as  a 
design  criterion 

Code  to  determine  type  of  wall  design, 
=  STD,  implies  standard  wall  sizes 

A,  B,  and  C  from  ASTM  C-76 
=  OPT,  implies  continuous  wall 

sizes  are  assumed  for  optimum 

design 
=  FIX,  implies  fixed  wall  size  is 

specified  and  only  steel  is 

designed 

Default  =  STD 

Cover  to  center  of  steel  reinforce- 
ment for  both  inner  and  outer  cages; 
or,  minimum  cover  of  elliptical  cage 
Default  =  1.25  in. 

Steel  area  ratio  of  outer-to-inner 
cages,  only  applies  to  RSHAPE  = 
CIRC, 

Default  =  0.75 


Note 
(7) 


(8) 


***   GO  TO  SECTION  C  *** 
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Reinforced  Concrete  Commentary 

(1)  For  round  pipe,  PDIA  is  the  pipe  diameter  measured  from  inner 
wall  to  inner  wall.  For  elliptical  pipe  (Level  2),  PDIA  is  the  average 
of  the  major  and  minor  diameters.   For  other  culvert  shapes  (Level 

3),  PDIA  is  a  nominal  average  diameter;  further  shape  definition  is 
required  in  Section  C. 

(2)  Specification  of  reinforcement  shape  by  "RSHAPE"  is  relative 
to  the  shape  of  the  concrete  wall.  That  is,  if  RSHAPE  =  CIRC,  two 
lines  of  reinforcement,  inner  and  outer,  run  parallel  to  the  inner 

and  outer  pipe  wall  periphery,  regardless  of  overall  pipe  shape.   If 
RSHAPE  »  ELLI,  one  line  of  reinforcement  traverses  the  pipe  wall  from 
an  interior  wall  position  at  crown  to  an  exterior  wall  position  at 
the  springline  and  back  to  an  interior  position  at  the  invert.  Lastly, 
if  RSHAPE  =  ARBI,  two  lines  of  reinforcement  can  be  specified  in  any 
shape  relative  to  the  pipe  wall.  Also,  the  steel  areas  may  be  selectively 
varied  around  the  pipe;  however,  this  option  only  applies  to  the  analysis 
mode. 

(3)  As  a  general  rule,  concrete  cracking,  nonlinear  compression, 
and  steel  yielding  should  be  used  for  all  problems  (NONLIN  ■  3) .  Lesser 
degrees  of  nonlinearity  are  useful  for  comparative  studies.  Note,  linear 
modeling  with  no  cracking  is  achieved  by  setting  NONLIN  =  1  and  STNMAT(1) 
■  large  value. 


25 


(4)  The  nonlinear  concrete  stress-strain  model  is  defined  with 
measures  of  strain;  i.e.,  cracking  strain  STNMAT(1),  compressive  strain 
at  elastic  limit  STNMAT(2),  and  strain  at  compressive  strength  STNMAT(3), 
Figure  II-3  illustrates  the  concrete  model  and  input  parameters. 

Reinforcing  steel  is  assumed  elastic  perfectly  plastic  and 
does  not  slip  relative  to  the  concrete  web  (i.e.,  deformed  bars  or 
mesh) . 

(5)  PFPC  (i.e.,  f')  and  other  parameters  on  this  card  are  supplied 
with  recommended  default  values  for  reinforced  concrete  pipe.  Thus, 

a  blank  data  card  may  suffice  for  this  input. 

(6)  Only  include  this  card  if  XMODE  »  ANALYS.  Definition  of 

the  parameters  ASI,  ASO,  TBI,  and  TBO  depend  on  RSHAPE.  For  the  special 
case,  RSHAPE  »  ARBI,  these  parameters  are  input  for  each  pipe  segment 
or  element.  Thus,  card  3B  is  repeated  ten  times  for  Level  1  or  2  and 
NPMAT  times  for  Level  3,  such  that  ASI  and  TBI  will  define  one  line 
of  reinforcement;  TBI  is  the  distance  from  the  inner  wall  to  the  center 
of  steel  area  ASI.  Similarly,  ASO  and  TBO  define  a  second  line  of  rein- 
forcement measured  from  the  outer  wall. 

For  the  more  common  cases,  RSHAPE  ■  CIRC  or  ELLI,  card  3B 
is  input  once  as  indicated. 

(7)  Only  include  this  card  if  XMODE  -  DESIGN.  Normal  ranges 
of  specified  safety  factors  are: 
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Steel  yielding,  1.5  -  2.0 
Concrete  crushing,  1.5  -  2.0 
Shear  failure,  2.0  -  3.0 
Crack  width,  0.01  in. 

(8)   In  the  design  mode  the  objective  is  to  determine  both  concrete 
wall  thickness  and  required  steel  area(s)  when  GOAL  ■  OPT  or  STD.  For 
the  case  GOAL  «  FIX,  the  specified  wall  thickness  is  held  constant, 
and  only  the  steel  area  is  designed.   In  this  case  only  the  crack  width 
and  steel  yielding  criterion  are  used  in  the  design  process. 


Tensile  rupture 


Figure  II-3.  Concrete  stress-strain  model. 
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Plastic  Pipe  Input  (Cards  IB,  2B) 

Card  IB.   Pipe  size  and  material  properties: 


Columns 

Variable 

(format) 

(units) 

01-10 

PDIA 

(F10.0) 

(in.) 

11-20 

PE 

(F10.0) 

(psi) 

21-30 

PNU 

(F10.0) 

31-40 

PULT 

(F10.0) 

(psi) 

41-50 

PDEN 

(F10.0) 

(pci) 

Entry  Description  Note 

Inside  pipe  diameter,  average  (1) 


Young's  modulus  of  pipe  material, 
Default  =  1,600,000  psi 

Poisson's  ratio  of  pipe  material, 
Default  =  0.35 

Ultimate  material  stress  at  rupture,     (2) 
Default  =  25,000  psi 

Density  of  pipe  material  used  for 
body  weight  loading;  only  applies  to 
Levels  2  and  3, 
Default  =0.0 
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Card  2B.   For  analysis  only,  plastic  pipe  wall  thickness: 


Columns 
(format) 

01-10 
(F10.0) 


Variable 
(units) 

PT 
(in.) 


Entry  Description 

Thickness  of  smooth  pipe  wall 


Card  2B.   For  design  only,  desired  safety  factors: 


Columns 
(format) 

01-10 
(F10.0) 


11-20 
(F10.0) 


21-30 
(F10.0) 


Variable 
(units) 

PFS(l) 


PFS(2) 


PFS(3) 


Entry  Description 

Desired  safety  factor  against 
excessive  deflection  (20%  of  PDIA) , 
Default  =4.0 

Desired  safety  factor  against  outer 
fiber  stress  rupture, 
Default  =3.0 

Desired  safety  factor  against 
elastic  buckling, 
Default  =3.0 


Note 
(3) 


Note 


(A) 


***  GO  TO  SECTION  C  *** 
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Plastic  Pipe  Commentary 

(1)  For  round  pipe,  PDIA  is  the  inside  diameter.  For  elliptical 
pipe  (Level  2),  PDIA  is  the  average  of  the  major  and  minor  inside  diameters. 
For  other  culvert  shapes  (Level  3),  PDIA  is  a  nominal  average  diameter; 
further  shape  definition  is  required  in  Section  C. 

(2)  The  stress-strain  law  for  plastic  pipe  is  assumed  linear 

up  to  brittle  failure  at  the  stress  level  PULT,  as  indicated  in  Figure 
II-4.  The  default  values  for  card  2B  are  representative  of  plastic 
reinforced  mortar  pipe.  For  other  types  of  brittle  plastic  the  user 
must  supply  the  appropriate  properties.  Note,  the  plastic  pipe  model 
is  suitable  for  other  brittle  materials,  such  as  cast  iron  or  clay 
pipe, 

(3)  Only  include  this  card  if  XMODE  -  ANALYS.  The  pipe  wall 
is  assumed  to  be  uniform  and  homogeneous  with  thickness  PT. 

(4)  Only  include  this  card  if  XMODE  -  DESIGN.  The  normal  range 
of  specified  safety  factors  are: 

Excessive  deflection,  3.0  -  4.0 
Outer  fiber  rupture,  2.5  -  3.5 
Elastic  buckling,  2.5  -  3.5 
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rupture 


Symmetric  in  tension 
or  compression 


Strain 


Figure  II-4.   Stress-strain  model  for  plastic. 
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Basic  Input,  (Cards  IB,  2B) 
Card  IB.   Pipe  size: 

Entry  Description 
Diameter  of  pipe,  average 

Card  2B.   Pipe  properties  (repeat  card  as  necessary) 


Columns 

Variable 

(format) 

(units) 

01-10 

PDIA 

(F10.0) 

(in.) 

Note 
(1) 


Columns 
(format) 

Variables 
(units) 

01-05 
(15) 

ISEQ1 

06-10 
(15) 

ISEQ2 

11-20 
(F10.0) 

PE 
(psi) 

21-30 
(F10.0) 

PNU 

31-40 
(F10.0) 

PA   2 
(in.  /in.) 

41-50 
(F10.0) 

PI  4 
(in.  /in.) 

51-60 
(F10.0) 

PDEN 
(lb/in.) 

Entry  Description 

First  pipe  element  material- 
sequence  number 

Second  pipe  element  material- 
sequence  number 

Young's  modulus  of  pipe  material 
for  the  sequence  (ISEQ1  to  ISEQ2) 

Poisson's  ratio  of  pipe  material 
for  the  sequence 

Area  of  pipe  section  per  unit  length 
for  the  sequence 

Moment  of  inertia  of  pipe  section  per 
unit  length  for  the  sequence 

Weight  of  pipe  element  per  unit  length 
for  the  sequence 


Note 
(2) 


(3) 


***  GO  TO  SECTION  C  *** 
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Basic  Commentary 

(1)  The  basic  pipe  model  allows  specification  of  pipe  properties 
that  are  nonuniform  around  the  periphery  of  the  culvert,  and,  thus, 

is  best  used  with  Levels  2  or  3.  Only  analysis  problems  can  be  considered, 
For  round  pipes,  PDIA  is  the  average  diameter  from  wall  center  to  wall 
center.  For  elliptical  shapes  (Level  2),  PDIA  is  the  average  of  major 
and  minor  diameters.  For  other  shapes  (Level  3),  PDIA  is  not  used; 
instead,  pipe  shape  is  defined  in  Section  C. 

(2)  For  Levels  2  and  3,  each  pipe  element  is  identified  with 

a  separate  material  number  in  a  consecutive  sequence  from  1  to  NPMAT, 
traveling  clockwise  around  the  pipe  (NPMAT  ■  number  of  pipe  elements, 
NPMAT  =  10,  for  Levels  1  and  2).   If  consecutive  groups  of  elements 
have  the  same  material  properties,  the  use  of  ISEQ1  and  ISEQ2  reduces 
the  number  of  2B  data  cards.  For  example,  if  all  pipe  elements  have 
the  same  properties,  then  only  one  card  is  required  with  ISEQ1  «  1 
and  ISEQ2  -  NPMAT. 

(3)  Linear  stress-strain  models  are  assumed  for  all  pipe  materials 
described  by  BASIC.  No  default  material  properties  are  provided,  and 
no  safety  factor  evaluation  is  determined. 
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SECTION  C  -  SOLUTION  LEVEL  DESCRIPTION 

Go  to  subsection  corresponding  to  chosen  solution  level: 

•  Level  1 

•  Level  2 

•  Level  3 

Commentaries  are  included  after  the  input  instructions  for  each 
solution  level.  Note,  Level  2  subsection  is  followed  by  an  additional 
subsection  for  "extended"  Level  2  operations.  Extended  Level  2  allows 
for  selective  modifications  of  standard  Level  2  mesh  models. 
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Level  1  Input  (Cards  1C,  2C) 

Card  1C.   Level  1  system  parameters 


Columns 

Variable 

(format) 

(units) 

Entry  Description 

Note 

01-10 
(F10.0) 

SDEN 
(pcf) 

Density  of  soil  over  pipe 

(1) 

11-15 
(15) 

NINC 

Total  number  of  load  steps  or  soil, 
lifts, 

Default  =  1 

(2) 

16-20 
(15) 

ISLIP 

Code  for  treatment  of  pipe-soil 
interface, 

(3) 

Card  2C, 


Columns 
(format) 

01-10 
(F10.0) 

11-20 
(F10.0) 


21-30 
(F10.0) 


=  0,  implies  no  slip,  rigid 
interface  bond 

=  1,  implies  full  slip,  friction- 
less  interface 

Level  1  soil  parameters  (repeat  card  for  each  load  step,  1=1 
to  NINC): 


Variables 
(units) 

HT(I) 
(ft) 

SEH(I) 
(psi) 


SVH(I) 


Entry  Description  Note 

Soil  height  above  crown  at  load 
step  I  (current  total  height) 

Young's  modulus  of  soil  in  vicinity     (4) 

of  pipe  at  load  step  I  (secant 

values) 

Poisson's  ratio  of  soil  in  vicinity 
of  pipe  at  load  step  I 


***  END  OF  INPUT  FOR  ONE  PROBLEM  *** 
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Level  1  Commentary 

(1)  Level  1  provides  a  solution  for  a  circular  pipe  deeply  buried 
in  a  homogeneous  soil.  All  loading  is  approximated  by  an  equivalent 
overburden  pressure  p  as  indicated  in  Figure  II-5.  To  simulate  gravity 
loads,  Level  1  uses  p  =  SDEN  [HT(I)  -  HT(I-1)]  for  each  load  step.  If 
it  is  desired  to  simulate  concentrated  loads,  the  user  must  artificially 
increase  the  soil  height  HT  by  an  amount  AH  to  represent  the  additional 
pressure  at  the  pipe  crown.  Pressure  at  the  crown  due  to  a  surface 
concentrated  load  Q  may  be  estimated  by  assuming  the  load  spreads  in 

a  30-degree  cone.  Therefore,  AH  »  3Q/[  (tt)  (SDEN)  (HT2)  ] .  For  additional 
details  see  Chapter  9  of  the  engineering  manual  [1], 

In  general,  gravity  loads  overshadow  contributions  from  concentrated 
loads  when  fill  heights  exceed  5  feet.  When  concentrated  loads  are 
significant  and  the  cover  depth  is  shallow,  Level  2  should  be  used. 

(2)  In  most  design  cases  it  is  adequate  to  place  the  entire  load 
in  one  step  (NINC  =  1).  However,  if  it  is  desired  to  study  the  response 
history  of  the  pipe,  five  to  ten  steps  should  be  used. 

(3)  Two  pipe-soil  interface  conditions  can  be  considered  -  slip 

or  no  slip.  No  slip  implies  both  normal  and  shear  forces  are  transmitted 
across  the  interface,  whereas  slip  only  transmits  normal  forces.  Generally, 
no  slip  is  recommended  for  design  problems. 
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(4)  As  a  rule,  the  Young's  modulus  of  soil  Increases  for  Increased 
states  of  confining  compression.  Accordingly,  specified  values,  SEH(I) , 
should  correspond  to  the  accumulated  overburden  pressure,  p  »  SDEN*HT(I) . 
Whenever  possible  specified  values  should  be  based  on  experimental 
tests  of  the  site  soil  or  other  direct  sources  of  information.  However, 
if  no  information  is  available,  Table  II-1  can  be  used  as  a  guide  to 
estimate  SEH(I)  as  a  function  of  overburden  pressure.  The  table  includes 
three  broad  soil  catagories  -  granular,  mixed,  and  cohesive  with  compaction 
ranging  from  fair  to  good.  Also  shown  are  estimates  for  Poisson's 
ratio  and  soil  density. 

If  several  load  steps  are  specified  (NINC  >  1),  Card  2C  is  repeated 
for  each  step,  wherein  SEH(I)  is  the  secant  modulus  for  accumulated 
fill  height,  HT(1).  Level  1  automatically  converts  the  secant  values 
into  incremental  chord  values  applicable  to  the  load  step.   See  Chapter 
6  of  the  engineering  manual  for  more  discussion  on  soil  models. 
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Overburden  Pressure,  Pf 


c 
c 
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D 
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D 
13 

D 
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I        I       I        y       y        f 


Homogeneous  Soil 


Es  =  Young's  Modulus 
vs  =  Poisson's  Ratio 


pipe 


A     ~T~  ~T~  ~T~  ~T~   ~1       F~  ~1 

Figure  II-5.   Approximate  model  representation  of  elasticity  boundary  value  problem. 
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Level  2  Input  (Cards  1C,  2C,  3C,  ID,  2D) 

Card  1C.   Define  mesh  type,  title,  and  special  options: 

Entry  Description 

Name  to  identify  type  of  automatic 

mesh, 

=  EMBA,  embankment  mesh 

=  TREN,  trench  mesh 

=  HOMO,  homogeneous  mesh 


Columns 

Variable 

(format) 

(units) 

01-04 

WORD 

(A4) 

05-72 
(17A4) 

73-76 
(A4) 


77-80 
(A4) 


TITLE 


W0RD1 


WORD  2 


User  description  of  mesh  to  be 
printed  with  output 

Command  to  include  frictional 
pipe-soil  interface  model, 

=  SLIP,  include  frictional  inter- 
face 

^  SLIP,,  do  not  include  interface 

Command  to  permit  user  to  selectively 

modify  the  automatic  mesh, 

=  MOD,  mesh  will  be  modified 
4   MOD,  mesh  will  not  be  modified 
(left  justified) 


Note 
(1) 


(2) 


(3) 
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Card  2C.   Define  print  options  and  mesh  parameters: 


Columns 
(format) 

01-05 
(15) 


06-10 
(15) 


11-15 
(15) 


Variable 
(units) 

IPLOT 


IWRT 


MGENPR 


Entry  Description 

Signal  to  create  a  plot  data  tape 
on  unit  10, 

=  0,  no  data  tape  created 

=  1,  create  data  tape 

Signal  to  print  out  soil  responses 

for  all  elements, 

=  0,  no  soil  response  printed  out 
=  1,  print  out  soil  responses 

Code  to  control  amount  of  print  out 
of  mesh  data, 

=  1,  minimal  printout;  just 

control  data 
=  2,  above,  plus  node  and  element 

input  data 
=  3,  above,  plus  generated  mesh 

data 
=  4,  maximal  printout  of  input  data 
Default  =  3 


Note 


16-20 
(15) 


NINC 


Number  of  construction  increments, 
=  -1,  combine  all  lifts  into 

one  monolith 
=  0,  used  for  data  check  only; 

all  data  is  read  but  not  executed 
=  N,  number  of  construction 

increments  to  be  executed, 

N  =  1  to  10 


(4) 


21-30 
(F10.0) 


31-40 
(F10.0) 


41-50 
(F10.0) 


RDIA 


HTCOVR 
(ft) 


DENSTY 
(pcf) 


Ratio  of  horizontal  to  vertical  pipe 
diameter  to  define  elliptical  pipe, 
Default  =1.0  (circle) 

If  WORD  =  EMBA  or  HOMO,  HTCOVR  is        (5) 
total  fill  soil  height  above  spring- 
line;  if  WORD  =  TREN,  HTCOVR  is  fill 
soil  height  above  top  of  trench 

Density  of  soil  above  truncated  mesh     (6) 
to  be  used  as  equivalent  overburden 
pressure 

continued 
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Card  2C.  continued 

Columns  Variable 

(format)  (units)      Entry  Description                     Note 

51-60  PRESUR       Uniform  surface  pressure  applied  on     (7) 

(F10.0)  (psi)        last  construction  increment 
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Card  3C.   Define  backpacking  or  trench  geometry: 


Columns 
(format) 

01-05 
(15) 


21-30 
(F10.0) 


31-40 
(F10.0) 


Variable 
(units) 

NUMPEL 


11-20       TPAC 
(F10.0)      (in.) 


TRNDEP 
(ft) 


TRNWID 
(ft) 


Entry  Description  Note 

Code  to  determine  amount  of  back-       (8) 

packing  for  embankment  mesh  only, 
=  0,  no  backpacking  material; 
this  zone  is  replaced  with 
fill  soil  and  bedding 
=  N,  the  number  of  consecutive 
elements  starting  at  the  crown 
to  be  assigned  backpacking 
properties,  remaining  elements 
adjacent  to  pipe  are  assigned 
to  fill  soil  zone  and  bedding, 
N  =  1  to  10 

Thickness  of  backpacking  zone 
around  the  pipe  (D/16  <TPAC<3D/16, 
where  D  is  pipe  diameter);  EMBA 
only, 

Default  =  D/12 

Depth  of  trench  measured  from  the       (9) 
pipe  invert  (TRNDEP>D/4,  where  D 
is  the  vertical  pipe  diameter); 
only  required  for  WORD  =  TREN 

Width  of  trench  (1.25D  <  TRNWID  < 
1.5D,  where  D  is  horizontal  pipe 
diameter);  only  required  for 
WORD  =  TREN 


Cards  4C-7C.   Cards  for  modifying  Level  2  meshes  go  here: 

If  W0RD2  =  MOD,  go  to  Extended  Level  2  for  data  input  instructions; 
otherwise  skip  Cards  4C-7C,  and  complete  input  with  cards  ID  and  2D. 
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Card  ID.   Material  identifier  card  (repeat  cards  ID  and  2D  for  each 
material) : 


Columns 
(format) 

01-01 
(Al) 


Variable 
(units) 

LIMIT 


02-05 
(14) 


06-10 
(15) 


ITYP 


11-20 
(F10.0) 


DEN  (I) 
(pcf) 


Entry  Description  Note 

Last  material  card-set  indicator; 
=  0,  read  another  set  of  material 

definitions 
=  L,  this  is  the  last  material 

input 

Material  identification  number;         (10) 
corresponds  to  the  material  zone 
of  the  chosen  mesh  type, 

=  1,  for  in-situ  soil  zones  or 

homogeneous  mesh 
=  2,  for  bedding  zones 
=  3,  for  embankment  fill  or  trench 

fill 
=  4,  for  embankment  backpacking 
or  trench  overfill 
(The  above  is  predefined  for  Level 
2  only . ) 

Selection  of  material  model  to  be       (11) 
associated  with  material  zone  I, 
=  1,  linear  elastic  (isotropic) 
=  2,  linear  elastic  (orthotropic) 
=  3,  (unspecified,  to  be  supplied 

by  user) 
=  4,  overburden  dependent 
=  5,  fully  nonlinear,  Extended- 

Hardin  model 
=  6,.  frictional  interface  model 
(If  WORD  1  =  SLIP,  see  comment) 

Density  of  material  I  used  to  compute 
gravity  loads;  only  applies  to 
ITYP  =  1,  2,  3,  4,  and  5. 


21-40 
(5A4) 


MATNAM 


User  description  of  material;  if 
ITYP  -  1,  2,  3,  4,  or  6,  MATNAM  may 
be  any  descriptive  information  to  be 
printed  out  with  data,  for  ITYP  =  5, 
MATNAM  is  used  to  select  the  soil 
classification  for  the  Extended- 
Hardin  model  as  follows  (left 
justified) : 


continued 
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Card  ID   continued 


Columns 

Variable 

(format) 

(units) 

21-40 

MATNAM 

(5A4) 

Entry  Description 

=  GRAN,  for  modeling  granular  soils 
=  MIXE,  for  modeling  mixed  soils 
=  COHE,  for  modeling  cohesive  soils 
=  TRIA,  for  utilizing  data  from  tri- 
axial  tests 
Default  =  MIXE 


Note 


***  GO  TO  CARD  2D  CORRESPONDING  TO  ITYP  *** 
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Card  2D.   ITYP  =  1,  linear  elastic: 


Columns 
(format) 

01-10 
(F10.0) 

11-20 
(F10.0) 


Variable 
(units) 

E 
(psi) 

GNU 


Entry  Description 

Young's  modulus  of  material  I 

Poisson's  ratio  of  material  I 


Card  2D.   ITYP  =  2,  orthotropic,  linear  elastic: 


Columns 
(format) 

01-10 
(F10.0) 

11-20 
(F10.0) 

21-30 
(F10.0) 

31-40 

(F10.0) 

41-50 
(F10.0) 


Variable 
(units) 

CP(1,1) 
(psi) 

CPU,  2) 
(psi) 

CP(2,2) 
(psi) 

CP(3,3) 
(psi) 

THETA 
(deg) 


Entry  Description 

Constitutive  parameter  at  matrix 
position  (1,1) 

Constitutive  parameter  at  matrix 
position  (1,2) 

Constitutive  parameter  at  matrix 
position  (2,2) 

Constitutive  parameter  at  matrix 
position  (3-3) 

Angle  of  the  material  axis  with 
respect  to  the  global  x-axis 


Note 
(12) 


Note 
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Card  2D.   ITYP  =  3,  user  developed  soil  model: 

Provisions  have  been  made  for  users  to  substitute  their  own  soil 
model  in  CANDE,  in  which  case  data  are  supplied  at  this  juncture.   See 
system  manual. 

Card  2D.   ITYP  =  4,  overburden  dependent  model,  user  table,  (repeat  card 

as  needed;  last  card  is  blank): 


Columns 
(format) 

Variable 
(units) 

01-10 
(F10.0) 

H(N) 
(psi) 

11-20 
(F10.0) 

E(N) 
(psi) 

21-30 
(F10.0) 

GNU(N) 

Entry  Description  Note 

Overburden  pressure  for  table  entry  N    (13) 


Young's  secant  modulus  for  table 
entry  N 

Poisson's  ratio  for  table  entry  N 


Ul 


Card  2D.   ITYP  =  5,  and  MATNAM  =  GRAN,  MIXE,  or  COHE;  Extended-Hardin 
model  for  three  soil  classifications: 


Columns 
(format) 

01-10 
(F10.0) 

11-20 
(F10.0) 

21-30 
(F10.0) 


Variable 
(units) 

XNUMIN 


XNUMAX 


XQ 


Entry  Description 

Poisson's  ratio  at  low  shear  strain, 
Default  =  0.10 

Poisson's  ratio  at  high  shear  strain, 
Default  =  0.49 

Shape  parameter  q  for  Poisson's  ratio 
function, 

Default  =  0.26 


Note 
(14) 


31-40 
(F10.0) 


VOIDR 


Void  ratio  of  soil,  range  0.1  to  3.0; 
see  comments  for  default 


41-50 
(F10.0) 

51-60 
(F10.0) 


SAT 


PI 


Ratio  of  saturation,  range  0.0  to  1.0; 
see  comments  for  default 

Plasticity- index/100,  range  0.0  to 
1.0;  see  comments  for  defaults 


61-65 
(15) 


NON 


Maximum  iterations  per  load  step, 
Default  =  5 


Card  2D.   ITYP  =  5,  and  MATNAM  =  TRIA;  Extended-Hardin  model  for 
triaxial  data  input 


Columns 
(format) 


Variable 
(units) 


Entry  Description 


Note 


01-30 
(3F10.0) 

31-40 
(F10.0) 


SI 


Same  as  card  above  (XNUMIN,  XNUMAX, 
XQ) 

Hardin  parameter  used  to  calculate 
maximum  shear  modulus 


41-50 
(F10.0) 


CI 


Hardin  parameter  used  to  calculate 
reference 


51-60 
(F10.0) 

61-65 
(15) 


NON 


Hardin  parameter  used  to  calculate 
hyperbolic  shear  strain 

Maximum  iterations  per  load  step, 
Default  =  5 
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Card  2D  ITYP  =  6,  interface  property  definition  (must  be  defined  with 
Card  ID  ten  times  if  WORD  1  =  SLIP): 


Columns 

Variable 

(format) 

(units) 

01-10 

ANGLE 

(F10.0) 

(deg) 

Entry  Description 

Angle  from  x-axis  to  normal  of 
interface  at  the  nodel  point 
(direction  of  normal  is  IX(1)  to 
IX(2);  see  comment) 


Note 
(15) 


11-20 
(F10.0) 


FCOEF 


Coefficient  of  friction  between  pipe 
and  soil  at  interface  nodal  point 


21-30 
(F10.0) 


TENSIL 
(lb) 


Tensile  breaking  force  at  contact  point, 
Default  =0.0 


***  INPUT  COMPLETED  FOR  ONE  PROBLEM  *** 
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Level  2  Commentary 

(1)  Level  2  allows  the  selection  of  three  fundamental  culvert 
installation  models;  the  embankment  mesh  (Figure  II-6) ,  the  trench 
mesh  (Figure  II-7) ,  and  the  homogeneous  mesh,  where  all  soil  zones 
have  the  same  material  properties;  the  construction  increments  are  as 

in  the  embankment  mesh.  Level  2  is  restricted  to  installations  symmetric 
about  the  vertical  centerline  of  the  pipe.  Accordingly,  only  half  the 
system  is  modeled  with  finite  elements.  The  automatic  generation  defines 
all  nodes,  elements,  construction  increments,  material  zones,  and 
boundary  conditions.  The  basic  mesh  topology  contains  86  quadrilateral 
elements  for  the  soil  zones  and  ten  pipe  elements  (see  Extended  Level 
2  for  mesh  topology) .  By  selectively  identifying  each  element  with 
a  material  zone  number  and  a  construction  sequence  number,  the  embankment 
mesh  and  trench  mesh  are  automatically  constructed  from  the  same  mesh 
pattern. 

(2)  If  W0RD1  «  SLIP,  the  mesh  is  automatically  altered  to  include 
eleven  interface  elements  at  common  nodes  between  pipe  and  soil.  The 
interface  model  allows  for  frictional  sliding,  complete  separation, 
and  rebonding  of  the  pipe-soil  interface  during  the  loading  schedule. 
This  option  should  be  restricted  to  special  problems. 

(3)  If  it  is  desired  to  alter  the  standard  Level  2  configurations, 
set  W0RD2  =  MOD.  In  so  doing,  the  user  must  supply  additional  data 
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J 


j  Surface  Loads 


I ][ :' 


P 


Fill  Soil  (3,5) 


■3 


Construction 
Lifts 


Zone  Description 
(a,  b) 

a  =  material  no. 
b  =  increment  no. 


■12  R, 


Figure  II-6.  Embankment  model  of  Level  2  boundary  value  problem. 
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Surface  loads 


Overfill 
Variable 


c 


c 


Trench 

Depth 

Variable 


3  Ri 


Overfill  Soil 
(4,5) 


3 
3 


<A    l     Construction 
'  Lifts 


Zone  description 
(a,b) 

a  =  material  no. 
b  =  increment  no. 


Figure  II-7.  Trench  model  of  Level  2  boundary  value  problem. 
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cards  C4  through  C7  when  indicated.  Some  motivations  for  altering  the 
standard  Level  2  mesh  are:   applying  concentrated  loads,  changing  bedding 
shape,  incorporating  an  imperfect  trench,  modifying  culvert  shapes, 
etc.  It  is  wise  to  exercise  standard  Level  2  options  before  using  the 
Extended  Level  2  option. 

(4)  A  maximum  number  of  ten  construction  increments  may  be  used 
to  simulate  the  construction  process  of  placing  soil  in  layers.  The 
first  increment  is  the  in-situ  soil,  bedding,  and  pipe;  second,  fill 
soil  up  to  the  springline;  third,  fill  from  springline  to  crown;  fourth, 
fill  from  crown  to  1.25  diameter  above  springline;  and  fifth,  fill 

up  to  a  maximum  of  2  diameters  above  springline.  Thereafter,  equivalent 
overburden  pressure  is  used  if  required. 

(5)  For  the  embankment  or  homogeneous  mesh,  HTCOVR  (see  Figure 
II-6)  specifies  the  total  fill  height  above  the  springline,  which  must 
be  at  least  0.75  diameter.   If  HTCOVR  is  less  than  2  diameters,  the 
mesh  surface  coincides  with  HTCOVR.  On  the  other  hand,  if  HTCOVR  is 
greater  than  2  diameters,  the  mesh  surface  will  be  truncated  at  2  diameters 
above  the  springline,  and  the  remaining  soil  weight  will  be  automatically 
applied  to  the  mesh  surface  in  equal  increments  of  overburden  pressure 

for  the  remaining  construction  increments. 

For  the  trench  mesh  (see  Figure  II-7) ,  HTCOVR  specifies  the  additional 
height  of  overfill  above  the  top  of  the  trench,  if  any.  If  the  combined 
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trench  fill  and  overfill  height  is  above  2  diameters  from  the  springline, 
the  mesh  will  be  truncated  at  2  diameters,  and,  the  remaining  soil  weight 
will  be  applied  in  equal  increments  of  overburden  pressure. 

(6)  When  the  mesh  surface  is  truncated  at  2  diameters  above  the 
springline,  the  user  must  specify  the  soil  density  of  the  remaining 
soil  above  the  mesh  surface  to  be  used  in  computing  the  additional 
overburden  pressure;  i.e.,  p  ■  density  (soil  height  -  mesh  height). 

(7)  To  mimic  Level  1  or  to  approximate  deep  fill  situations,  surface 
pressure  loads,  PRESUR,  can  be  used  in  place  of  gravity  loads.  In  such 
cases,  set  NINC  ■  -1,  DENSTY  -0.0,  and  determine  PRESUR  as  the  product 
of  fill  height  above  crown  times  soil  density. 

(8)  The  automated  backpacking  feature  is  only  applicable  to  the 
embankment  mesh  (see  Figure  II-6).  (For  the  trench  mesh,  Extended 
Level  2  can  be  used  to  define  backpacking.)  There  are  ten  candidate 
backpacking  elements  immediately  adjacent  to  the  ten  pipe  elements. 
The  variable  NUMPEL  specifies  how  many  consecutive  elements  starting 
at  the  crown  are  assigned  backpacking  properties.  For  example,  if 
NUMPEL  =  5,  backpacking  covers  the  top  half  of  the  pipe;  if  NUMPEL 

«■  10,  the  entire  pipe  is  surrounded  with  a  backpacking  ring,  including 
part  of  the  bedding  zone. 
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(9)  The  trench  depth   (TRNDEP)   is  automatically  scaled  up  to  the 
nearest  quarter  diameter  depth.   For  example,  if  the  trench  depth  is 
specified  between  1.0  and  1.25  diameters,   the  trench  depth  is  set  at 
1.25   diameters.  For  trench  depths   greater  than   1.25   diameters,   the 
mesh  surface  coincides  with  the  top   of  the   trench,   and  any  overfill 
soil   (material  zone  4)   is  treated  as  equivalent  overburden  pressure. 
To  avoid   **skinny"  elements,  do  not  use  trench  depths  in  the   range 
1.25  to  1.30  diameters. 

(10)  The  user  is   free  to  assign  any  material  stress -strain  model, 
ITYP  -   1,  2,   3,   A,  or  5,   to  each  material  zone   I.   The  material  zone 
numbers  correspond  to  predefined  areas  of  the  Level  2  mesh.   For  the 
homogeneous  mesh    (WORD  «  HOMO) ,   the  entire  mesh  is   zoned  with  I  =   1 . 
For  the  embankment  mesh,   the  zones   are:   in-situ  soil   (I  =   1),  bedding 
(I  -  2) ,   fill  soil   (I  -  3) ,   and  backpacking    (I  -  4) .   Backpacking  only 
needs  to  be  specified  if  NUMPEL  »  0.   Lastly,   for  the  trench  mesh,   the 
zones  are:   in-situ  soil   (I  ■  1),  bedding   (1-2),  trench  fill    (I  = 
3),  and  trench  overfill   (1=4). 

(11)  The  selection  of  the  material  stress -strain  model,   ITYP, 
depends  on  the  desired  characteristics  of  material  behavior.   A  full 
discussion  of  each  model  is  given  in  the  engineering  manual.   Naturally, 
the  model  selected  should  be  commensurate  with  experimental  information 
of  the  material  behavior.    If  information  is   scanty,   the  linear  model, 
ITYP  ■  1,  will  provide  reasonable  results.  A  more  realistic  model  for 
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soils  is  provided  with  overburden  dependent  model,  ITYP  =  4;  it  is 
recommended  for  general  use.  When  a  complete  set  of  triaxial  soil  data 
is  available,  the  fully  nonlinear  model,  ITYP  »  5,  provides  the  most 
realistic  soil  model.  Provision  has  been  made  for  the  user  to  incorporate 
his  own  soil  model  subroutine  in  CANDE  under  the  heading  ITYP  »  3. 

Lastly,  the  interface  model,  ITYP  »  6,  is  only  associated  with 
interface  elements  and,  for  Level  2  operations,  needs  only  to  be  input 
when  W0RD1  »  SLIP.  In  such  cases,  every  Interface  element  with  a  different 
normal  angle  is  identified  with  a  different  material  indicator,  I  (see 
comment  15). 

(12)  The  linear  elastic  model  implies  the  soil  zone  is  assigned 
constant  values  for  Young *s  modulus  and  Poisson's  ratio.  See  Table 
II- 1  for  typical  soil  properties. 

(13)  For  the  overburden  dependent  model,  Card  2D  is  repeated  two  to 
ten  times,  thereby  constructing  a  table  of  material  properties  versus 
overburden  pressure  for  material  zone  I.  (The  last  card  is  blank  to 
terminate  reading.)  Once  the  table  is  established,  CANDE  computes  the 
overburden  pressure  for  each  element  in  zone  I  by  adding  the  vertical 
stress  from  the  previous  load  step  to  one-half  of  the  current  increment 
of  body  weight  stress  above  each  element.  Using  this  value  of  overburden 
pressure,  the  material  properties  are  linearly  interpolated  from  the 
input  table  to  be  used  in  the  current  solution  step.  Material  properties 
are  assumed  constant  beyond  the  end  points  of  the  table. 
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Typical  overburden  dependent  soil  properties  are  given  in  Table 
II-1  and  can  be  used  if  no  other  data  are  available. 

(14)  The  Extended -Hardin  soil  model  employs  a  variable  shear  modulus 
and  a  variable  Poisson's  ratio.  The  user  should  consult  Chapter  6  of 
the  engineering  manual  before  implementing  this  model.  The  model  is 
best  used  in  conjunction  with  a  set  of  triaxial  data  (MATNAM  »  TRIA)  ; 
however,  the  model  can  be  defined  without  triaxial  tests  by  specifying 
the  soil  classification  (MATNAM  =  GRAN,  MIXE,  or  COHE),  along  with 
the  void  ratio,  percent  saturation,  and  plasticity  index. 

Default  values  for  these  quantities  are: 

GRAN      MIXE      COHE 
Void  ratio  0.60 

Percent  saturation/ 100   0.00 
Plasticity  index/ 100     0.00 

(15)  In  the  case  of  Level  2  with  W0RD1  =  SLIP,  the  material  properties 
for  the  interface  elements  are  different  for  each  element  because 

the  interface  normal,  ANGLE,  changes  around  the  pipe  periphery.  Therefore, 
Cards  1D  and  2D  must  be  input  for  each  of  the  eleven  interface  elements, 
starting  with  I  =  1  at  the  pipe  invert  counterclockwise  to  I  =  11  at 
the  crown.  If  the  pipe  is  circular,  the  respective  ANGLE  for  the  element 
normals  are:  -90,  -72,  -54,  -36,  -18,  0,  18,  36,  54,  72,  and  90  degrees. 


0.50 

1.00 

0.50 

0.90 

0.05 

0.20 
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If  the  interface  shear  force  exceeds  the  product  of  the  normal 
force  and  the  friction  coefficient  (FCOEF) ,  the  interface  element 
permits  relative  movement  according  to  the  Coulomb  friction  hypothesis. 
Should  the  interface  normal  force  exceed  the  tensile  breaking  limit 
(TENSIL) ,  the  contact  surfaces  will  separate  from  each  other  and  only 
rebond  if  subsequent  loading  brings  the  interface  back  together.  Note, 
in  Level  3  operations,  interfaces  between  two  soil  bodies  can  be  defined 
as  well  as  interfaces  between  pipe  and  soil. 
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Extended 

Card  AC. 

Columns 
(format) 

01-05 
(15) 

06-10 
(15) 

11-15 
(15) 


Level  2  Input  (Cards  4C,  5C,  6C,  and  7C) 

Number  of  quantities  to  be  modified: 


Variable 
(units) 

NEWXY 


NEWEL 


NEWBD 


Entry  Description  Note 

Number  of  nodal  points  to  be  changed     (1) 
with  card  5C 

Number  of  element  definition  cards  to 
be  changed  with  Card  6C 

Number  of  new  boundary  conditions 
that  will  be  defined  with  Card  7C, 
note,  this  number  is  restricted  to 
less  than  50  due  to  storage  require- 
ments 


Redefine  nodal  coordinates.   (Repeat  this  card  NEWXY  times; 


Card  5C. 

Redefine  n< 

skip  this 

Columns 

Variable 

(format) 

(units) 

01-05 

NEWNP 

(15) 

11-20 

XC00RD 

(F10.0) 

(in.) 

21-30 

YCOORD 

(F10.0) 

(in.) 

Entry  Description 

Number  of  nodal  points  whose 
coordinates  are  to  be  changed 

New  x-coordinate  location  of  NEWNP 


New  y-coordinate  location  of  NEWNP 


Note 
(2) 
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Redefines  element  parameters  (repeat  NEWEL  times;  skip  if 


Card  6C. 

Redefines 

NEWEL  =  0) 

Columns 

Variable 

(format) 

(units) 

1-5 

NE 

(15) 

6-10 

NP(1) 

(15) 

11-15 

NP(2) 

(15) 

16-20 

NP(3) 

(15) 

21-25 

NP(4) 

(15) 

26-30 

NP(5) 

(15) 

31-35 

NP(6) 

(15) 

36-40 

NP(7) 

(15) 

Entry  Description 

Element  number  whose  parameters 
are  to  be  changed 

New  node  number  I  of  element 
connectivity;  leave  blank  if  no 
change 

New  node  number  J  of  element 
connectivity;  leave  blank  if  no 
change 

New  node  number  K  of  element 
connectivity;  leave  blank  if  no 
change 

New  node  number  L  of  element 
connectivity;  leave  blank  if  no 
change 

New  material  identification  number; 
leave  blank  if  no  change 

New  construction  increment  number; 
leave  blank  if  no  change 

Leave  this  entry  blank  unless  the 
interface  elements  are  to  be  altered 


Note 
(3) 
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Card  7C.   Define  additional  boundary  conditions  (Repeat  this  card  NEWBD 
times;  if  NEWBD  =  0,  skip  this  card): 


Columns 
(format) 

01-05 
(15) 

06-10 
(15) 


Variable 
(units) 

NU 


IFLAG(l) 


Entry  Description 

Nodal  point  number  where  additional 
boundary  condition  is  specified 

Boundary  condition  code  for 
x-coordinate: 

=  0,  force  boundary  condition 

(load) 
=  1,  displacement  boundary 
condition 


Note 
(4) 


11-20 
(F10.0) 


21-25 
(15) 


BV(1)        Applied  x-force  load  (lb/in.)  or 
(lb/in.      x-displacement  (in.),  depending  on 
or  in.)      value  of  IFLAG(l) 

IFLAG(2)     Boundary  condition  code  for 
y-coordinate: 

=  0,  force  boundary  condition 

(load) 
=  1,  displacement  boundary 
condition 


26-35 
(F10.0) 


36-45 
(F10.0) 


46-50 
(15) 


BV(2) 
(lb/in. 
or  in. ) 

BY  (3) 
(deg) 


IA 


Applied  y-force  load  (lb/in.),  or 
y-displacement  (in.),  depending  on 
value  of  IFLAG(2) 

Angle  for  skew  boundary  conditions; 
angle  is  measured  counterclockwise 
from  the  global  x-axis 

Construction  increment  in  which 
boundary  condition  is  activated; 
force  boundary  conditions  are 
applied  only  during  increment  IA, 
Default  =1.0 


***  RETURN  TO  SECTION  C  FOR  CARDS  ID  AND  2D  *** 
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Extended  Level  2  Commentary 

(1)  Extended  Level  2  permits  the  user  to  adjust  the  standard  Level 
2  finite  element  mesh  data  to  redefine  nodal  coordinates,  to  add  or 
alter  material  zones,  to  change  the  construction  sequence  numbering, 
and  to  add  concentrated  loads  or  other  boundary  conditions. 

To  utilize  Extended  Level  2  effectively,  it  is  necessary  to  know 
the  mesh  configurations.  The  nodal  point  and  element  numbering  scheme 
is  the  same  for  all  meshes,  except  when  interface  elements  are  included 
(W0RD1  =  SLIP);  in  this  case  the  basic  numbering  scheme  is  altered 
to  accommodate  the  22  extra  nodes  required  by  the  11  interface  elements. 
Figures  II-8  and  II-9  illustrate  the  basic  element  and  node  numbering 
schemes.  The  alteration  of  the  basic  node  numbering  scheme  to  include 
interface  elements  is  identified  in  Table  II-2. 

For  either  configuration,  the  nodal  coordinates,  material  zones, 
construction  increment  numbering,  and  loading  depend  on  the  parameters 
specified  on  Cards  1C  through  3C.  Consequently,  the  user  should  obtain 
mesh  data  output  from  standard  Level  2  to  serve  as  reference  for 
Extended  Level  2. 

(2)  New  coordinates  can  be  specified  for  any  and  all  existing 
nodes;  see  Figure  II-9  for  reference.  Note  automatic  mesh  checking 
routines  are  by-passed  in  Extended  Level  2  operations;  therefore,  extreme 
care  must  be  taken  to  avoid  defining  nodal  coordinates  that  result 

in  "badly  shaped'*  or  ** inside-out "  elements. 
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Figure  II-8.   Basic  element  numbering  scheme. 
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Figure  11-9.    Basic  node  numbering  scheme. 
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Table  II-2.   Corresponding  Node' Table  for  W0RD1  =  SLIP 


Basic 

Slipping 

Basic 

Slipping 

Scheme 

Scheme 

Scheme 

Scheme 

Node 

Node(s)a 

Node 

Node(s)a 

1 

1 

50 

62 

2 

2 

51 

63 

• 

• 

• 

• 

• 

37 

• 

37 

• 

55 

• 

67 

38 

38,  39,  40 

56 

68,  69,  70 

39 

41,  42,  43 

57 

71,  72,  73 

40 

44,  45,  46 

58 

74,  75,  76 

41 

47,  48,  49 

59 

77,  78,  79 

42 

50,  51,  52 

60 

80,  81,  82 

43 

53 

61 

83 

44 

• 

54 

• 

62 

• 

84 

• 

• 

48 

• 

58 

• 
• 

• 

49 

59,  60,  61 

110 

132 

Instead  of  one  node  at  the  pipe-soil  interface  locations, 
three  nodes  are  defined  at  the  same  location  to  form  the 
interface  element.   For  every  node  triplet  above  (a,  b,  c) : 

a  =  pipe  node  IX(1) 

b  =  soil  node  IX(2) 

c  =  "free  node"  IX(3) 

Note:  The  eleven  interface  elements  are  numbered  97  through 
107,  beginning  at  the  invert  and  continuing  counterclockwise 
to  the  crown. 
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(3)  New  element  properties  can  be  defined  for  any  and  all  existing 
elements;  see  Figure  II-8  for  reference.  The  element  connectivity  array, 
NP(1),  NP(2),  NP(3),  and  NP(4) ,  contains  the  nodes  connected  to  each 
element  in  counterclockwise  order.  Generally,  there  will  be  little 
motivation  to  change  the  connectivity.  However,  the  material  identification 
number,  NP(5),  and  the  construction  increment  number,  NP(6) ,  allow 
selective  changes  to  material  zoning.  For  example,  bedding  zones  can 

be  expanded  by  assigning  those  elements  in  the  expanded  bedding  zone 
NP(5)  «=  2  and  NP(6)  =  1.  Further,  new  material  zones,  such  as  an  imperfect 
trench  zone  (e.g.,  NP(5)  =  5),  can  replace  existing  zones  as  desired. 
If  it  is  desired  to  set  any  NP(I)  value  to  zero,  enter  NP(I)  =  -1. 

(4)  The  standard  Level  2  options  are  limited  to  gravity  loads 
and  uniform  surface  pressure  loads.  A  prime  purpose  of  Extended  Level 
2  is  to  permit  consideration  of  concentrated  strip  loading,  which 

is  inherent  for  plane  strain,  by  assigning  to  it  any  node  (or  nodes) 
desired.  It  must  be  remembered  that  Level  2  employs  symmetry;  therefore, 
any  assigned  load  is  automatically  applied  symmetrically.  As  an  example, 
to  apply  a  concentrated  strip  load,  q  «■  1,000  lb/in. ,  on  the  soil  surface 
directly  over  the  pipe  crown,  set  NU  =  103  (basic  mesh)  and  BV(2)  ■  -500.0 
(i.e.,  -q/2);  all  other  entries  are  zero  (or  blank),  except  for  IA,  which 
defines  the  desired  load  increment  number. 

To  represent  a  single  concentrated  point  load,  Q,  by  a  strip  load, 
q,  let  q  -   3Q/4L,  where  L  is  the  shortest  distance  from  the  point  load 
to  the  pipe  crown.  This  approximation  is  based  on  equating  vertical 
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stresses  (at  the  pipe  crown)  between  free  field  theories  for  concentrated 
and  strip  loadings. 
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Level  3  Input  (Cards  1C,  2C,  3C,  4C,  5C,  ID,  2D) 
Card  1C.   Mesh  identification  card: 


Columns 

Variable 

(format) 

(units) 

01-04 

WORD 

(A4) 

(word) 

06-73 
(17A4) 


TITLE 
(words) 


Entry  Description 

Word  to  denote  user-defined  mesh 
is  forthcoming, 

=  PREP,  Level  3  data  are  read 

H  PREP,  error 

User  description  of  the  mesh  to  be 
printed  with  output 


Note 
(1) 
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Card  2C. 

Control  pa 

Columns 

Variable 

(format) 

(units) 

01-05 

NINC 

(15) 

06-10 

MGENPR 

(15) 

Control  parameters  for  Level  3  solution: 


Entry  Description 


Number  of  construction  increments; 
limits  of  NINC  are  1  to  10, 
Default  =  1 

Control  for  amount  of  printout  of  mesh 
data, 

=  1,  minimal  printout;  just  control 

data 
=  2,  above,  plus  node  and  element 

input  data 
=  3,  above,  plus  generated  input 

data 
=  4,  maximal  output  of  mesh  data 
Default  =  3 


Note 
(2) 


11-15 

NPUTCK 

(15) 

16-20 

IPLOT 

(15) 

21-25 
(15) 


26-30 
(15) 

31-35 
(15) 

36-40 
(15) 


IWRT 


NPT 


NELEM 


NBPTC 


Determines  if  execution  or  data  check: 
=  0,  execute  data 
=  1,  check  data  only,  terminate 

If  a  nonzero  value  is  entered,  the  plot- 
ting file,  unit  10,  will  contain  the 
input  mesh  data  and  structural  responses 
for  each  load  step 

If  a  nonzero  value  is  entered,  the 
interface  and  continuum  element 
responses  will  be  printed  out 

Total  number  of  nodal  points  in 
the  mesh 

Total  number  of  elements  in  the  mesh 


Number  of  boundary  conditions  for  this 
problem 
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Card  3C. 

Nodal  poin 

Columns 

Variable 

(format) 

(units) 

01-01 

LIMIT 

(Al) 

(letter) 

definition,  basic  generation  (repeat  as  needed) 


Entry  Description 

Signal  to  indicate  the  last  node 

input, 

\  L,  more  nodal  cards  to  come 
=-  L,  this  is  last  nodal  card 


Note 


02-05 
(14) 


NNP 


Nodal  point  number;  must  be  between 
1  and  NPT 


10-10 
(ID 


MODEG 


Defines  generation  mode  between 
previous  card  and  this  card, 

=  0,  standard  input,  no  generation 
=  2,  generates  nodes  between 

previous  card  and  current  card  in 
the  manner  requested  below 


(3) 


11-20 

XCOORD 

(F10.0) 

(in.) 

21-30 

YCOORD 

(F10.0) 

(in.) 

31-35 

NPINC 

(15) 

41-50 

SPACNG 

(F10.0) 

x-coordinate  of  NNP 


y-coordinate  of  NNP 


Increment  to  be  added  to  each  genera-    (4) 
ted  node,  if  MODEG  =  2  (positive), 
Default  =  1 

Ratio  of  distance  between  successively   (5) 
generated  pairs  of  nodes,  if  MODEG  =  2; 
only  applies  to  straight  line  genera- 
tion; 

Default  =1.0 


51-60 
(F10.0) 


RADIUS 
(in.) 


Defines  radius  of  generated  curve,  if   (6) 
MODEG  =  2, 

=  0.0,  straight  line  generation 
^  0.0,  circular  arc  generation; 
convexity  to  the  right  in  travel- 
ing from  previous  node  to  current 
node 
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Card  3C.   Nodal  point  definition,  advanced  generation  (repeat  as  needed) 

Column  9     Variable 

(format)     (units)      Entry  Description  Note 

01-01       LIMIT        Last  nodal  card  indicator 

(Al)        (letter)        *?  L,  implies  there  are  more  nodal 

cards 
=  L,  implies  this  is  the  last  nodal 

card 

02-05       NNP  Nodal  point  number  must  be  between  1 

(14)  and  NPT 

06-08       KRELAD       Permits  defining  the  coordinates  of     (7) 
(13)  the  node  NNP  by  referring  to  pre- 

viously defined  nodes  (associated 
with  MODEG  =  0,  or  2), 

=  0,  no  reference  to  another 

node 
=  1,  x-coordinate  will  be  specified 

by  a  previously  defined  node 
=  2,  y-coordinate  will  be  specified 

by  a  previously  defined  node 
=  3,  both  coordinates  will  be 
specified  by  previously  defined 
nodes 

09-09       LGTYPE       Allows  for  Laplace  generation  in  one     (8) 
(II)  coordinate  direction  only;  or  genera- 

tion along  a  quarter  of  an  ellipse, 
=  0,  no  action 
=  1,  x-coordinate  of  the  generated 

node  will  be  Laplace-generated 
=  2,  y-coordinate  of  the  generated 

node  will  be  Laplace-generated 
=  4,  one  quarter  of  an  ellipse  will 
be  generated 
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Columns     Variable 
(format)     (units) 


Entry  Description 


Note 


10-10       MODEG        Defines  generation  mode  between 
(II)  previous  card  and  this  card, 

=  0,  standard  input,  no  generation 

=  1,  no  generation;  recalls 
coordinates  of  a  previously 
defined  node 

=  2,  generates  nodes  between 
previous  card  and  this  card 

=  3,  generates  nodes  between  pre- 
vious card  and  this  card,  where 
the  nodal  point  NNP  on  this  card 
has  been  previously  defined 

-  5,  nonsequential  node  number 
input;  allows  for  the  genera- 
tion of  geometric  lines  with- 
out the  need  for  incremental 
node  numbering 

11-20       XCOORD       Defines  the  x-coordinate  of  NNP, 
(F10.0)      (in.  or         =  x-coordinate  value,  if  KRELAD 
node)  =0,2 

=  a  node  number  if  KRELAD  =  1  or  3 
(right  justified) 

21-30       YCOORD       Defines  the  y-coordinate  of  NNP, 
(F10.0)      (in.  or         =  y-coordinate  value,  if  KRELAD 
node)  =  0  or  1 

=  a  node  number,  if  KRELAD  =  2  or  3 
(right  justified) 


(9) 


31-35 
(15) 


NPINC 


Increment  to  be  added  to  each  generated 
node,  if  MODEG  =  2  or  3, 
Default  =1.0 


41-50 
(F10.0) 


SPACNG 


Ratio  of  distance  between  successively 
generated  pairs  of  nodes,  if  MODEG  =  2 
or  3;  only  applies  to  straight  line 
generation, 

Default  =1.0 


51-60 
(F10.0) 


RADIUS 


Defines- radius  of  generated  curve, 
=  0.0,  straight  line  generation 
V  0,  circular  arc  generation; 

convexity  to  the  right  in  travel- 
ing from  previous  node  to  current 
node 
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Card  4C.   Element  definition  card: 


Columns 
(format) 

01-01 
(Al) 


Variable 
(unit) 

LIMIT 
(letter) 


Entry  and  Description 

Last  element  card  indicator, 
^  L,  implies  there  are  more 

element  cards 
=  L,  implies  this  is  last 
element  card 


Note 


02-05 
(14) 

NE 

06-10 
(15) 

IX(1) 

11-15 
(15) 

IX(2) 

16-20 
(15) 

IX  (3) 

21-25 
(15) 


IX(4) 


Element  number;  must  be  between 
1  and  NELEM 

Node  number  I  of  element  connectivity 


Node  number  J  of  element  connectivity 


Node  number  K  of  element  connectivity, 
=  K,  node  for  quad  and  triangular 

elements 
=  0,  for  beam  elements 
=  A  node  unused  by  any  other  element 

for  interface  elements 

Node  number  L  of  element  connectivity, 
=  L,  node  for  quad  elements 
=  0,  for  beam,  triangular,  and 
interface  elements 


(10) 


26-30 
(15) 


IX  (5) 


Material  identification  number, 

=  1  to  10,  for  quad  and  triangular 

elements 
=  1  to  30,  for  beam  elements 
=  1  to  30,  for  interface  elements 


(11) 


31-35 
(15) 

IX(6) 

36-40 
(15) 

IX  (7) 

41-45 
(15) 

INTRAL 

Construction  increment  number  in  which 
element  enters  system, 
Default  =  1 

Set  to  1  to  designate  an  interface 
element  otherwise  leave  blank 

Node  increment  to  be  added  to  connect-   (12) 
ivity  nodes  for  element  row  generation, 
Default  =  1 


continued 
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Card   4C.      continued 


Column  s 
(format) 

Variable 
(unit  s) 

Entry  Description 

Note 

46-50 
(15) 

NUMLAY 

Number  of  element  rows  to  be 
generated, 

Default  =  1 

(13) 

51-55 
(15) 

INTERL 

Node  number  increment  between  element 
rows 
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Card  5C.   Boundary  conditions: 


Column  s 
(format) 

01-01 
(Al) 


Variable 
(units) 

LIMIT 


Entry  Description 

Last  boundary  condition  indicator, 

k   L,  implies  there  are  more  boundary 

condition  cards 
=  L,  implies  this  is  last  boundary 

condition  card 


Note 


02-05 
(14) 

06-10 
(15) 


NP  Nodal  point  number  for  boundary 

condition 

IIFLG(l)     Boundary  code  for  x-coordinate, 
=  0,  force  boundary  condition 
=  1,  displacement  boundary 
condition 


(14) 


11-20 
(F10.0) 


21-25 
(15) 


BIVD(l)      Value  of  boundary  condition  in  x 
direction, 

=  applied  force,  if  IIFLG(l)  =  0 
=  displacement,  if  IIFLG(l)  =  1 

IIFLG(2)     Boundary  code  for  y-coordinate 

=  0,  force  boundary  condition 
=  1,  displacement  boundary 
condition 


26-35 
(F10.0) 


36-45 
(F10.0) 

46-50 
(15) 


BIVD(2) 
(lb/in.) 


THETA 
(deg) 

IA 


Value  of  boundary  condition  in  y 
direction, 

=  applied  force,  if  IIFLG(2) 
=  0 

=  displacement,  if  IIFLG(2)  =  1 

Angle  for  skewed  boundary  condition; 
measured  counterclockwise  from  x-axis 

Construction  increment  in  which  this 
boundary  condition  is  activated;  force 
boundary  conditions  are  only  applied 
during  this  increment, 
Default  =  1 


(15) 


51-55 
(15) 


NNP 


Final  node  number  in  a  sequence  of 
boundary  conditions  to  be  generated 


(16) 


continued 


75 


Card  5C.   continued 


Columns     Variable 
(format)     (units) 


Entry  Description 


Note 


56-60 

INCR 

(15) 

61-70 

PJ 

(F10.0) 

(psi) 

71-80 

PK 

(F10.0) 

(psi) 

Node  number  increment  used  to 
generate  boundary  condition  nodes 
NP  to  NPP 

Pressure  magnitude  at  the  first  node 
of  the  sequence,  NP 

Pressure  magnitude  at  the  last  node 
in  the  sequence,  NNP 


(17) 
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Card  ID.   Material  identifier  card: 

Follow  instructions  for  Card  ID  of  Level  2  input  section.   The  only 
difference  is  that  the  material  number,  I,  corresponds  to  the  user's 
material  zone  scheme,  IX(5),  and  not  the  Level  2  material  zoning  scheme. 

Card  2D.   Material  parameters  (depends  on  ITYP) : 

Follow  instructions  for  Card  2D  of  Level  2  input  section  with  the 
corresponding  ITYP. 


***  INPUT  COMPLETED  FOR  ONE  PROBLEM  *** 
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Level  3  Comment  airy 

(1)  Level  3  is  the  customary  method  of  finite  element  input. 
Accordingly,  the  user  must  prepare  finite  element  mesh  data  representa- 
tive of  the  soil-culvert  system  to  be  designed  or  analyzed.  In  so  doing, 
four  groups  of  data  cards  must  be  completed  for  Level  3  in  addition 

to  the  Control  Cards  1C  and  2C.  These  cards  are:  (1)  nodal  point  definition 
(Cards  3C  basic  or  advanced),  (2)  element  definition  (Cards  AC),  (3) 
boundary  conditions  (Cards  5C) ,  and  (4)  material  property  definition 
(Cards  1D  and  2D,  same  as  Level  2). 

To  aid  the  user,  CANDE  is  equipped  with  many  sophisticated  mesh 
generation  features  that  can  greatly  reduce  the  amount  of  input  data. 
These  features  are  discussed  as  they  arise  in  the  input  instructions. 

(2)  The  number  of  construction  increments,  NINC,  dictates  how 
many  load  steps  CANDE  will  process.  All  elements  with  construction 
increment  numbers  less  than  or  equal  to  the  current  load  step  are  in- 
cluded in  the  current  stiffness  matrix. 

(3)  The  nodal  point  definition,  Card  3C,  is  given  with  two  versions 
to  provide  an  option  of  **basic"  or  ** advanced"  node  generation  schemes. 
Actually,  the  basic  option  is  a  subset  of  the  advanced  option,  but 

is  given  separately  to  minimize  confusion.  The  basic  option  provides 

node  generation  features  commonly  employed  in  most  finite  element  programs 
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and  is  relatively  straightforward.  On  the  other  hand,  the  advanced 
option  is  more  powerful  but  requires  some  ' 'learning' »  on  the  part 
of  the  user. 

For  the  basic  option,  when  MDDEG  ■  0,  the  input  node  NNP  is 
assigned  the  coordinates  XC00RD  and  YC00RD  and  other  entries  are  ig- 
nored. When  MODEG  =  2,  nodes  are  generated  beginning  from  the  node 
number  defined  on  the  previous  card  to  the  node  number  NNP  defined 
on  the  current  card. 

(4)  The  generated  nodes  are  numbered  by  successively  adding  NPINC 
to  the  previous  node  number  until  the  sum  is  greater  or  equal  to  NNP. 
The  last  node  is  always  numbered  NNP.  Accordingly,  the  number  of  nodal 
points  generated  is  (NNP  -  NNP*) /NPINC,  where  NNP*  is  the  node  number 
on  the  previous  card.  Node  generation  can  go  in  reverse,  i.e.,  NNP* 
can  be  greater  than  NNP,  but  NPINC  is  input  positive.  Nodes  can  be 
input  in  any  order. 

(5)  The  coordinates  assigned  to  generated  nodes  between  NNP*  and 
NNP  depend  on  SPACNG  for  straight  lines  and  RADIUS  for  circular  areas. 
If  SPACNG  -  1.0,  all  generated  nodes  will  be  evenly  spaced.  If  SPACNG 
>  1.0,  the  lengths  between  successively  generated  nodes  will  grow  with 
the  ratio  SPACNG.  Similarly,  if  SPACNG  <  1.0,  successive  lengths  will 
shrink  by  the  ratio  SPACNG. 
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(6)  To  generate  nodes  in  a  circular  arc,  specify  the  arc  radius. 
All  arc  segments   are  uniformly  spaced,   and  the   arc  generation  must 

be   less  or  equal  to   180  degrees.    If  RADIUS  is   specified  positive,   the 
convexity  of  the  arc  is  to  the  right  in  traveling  from  NNP*  to  NNP. 
Opposite  curvature  is   given  for  a  negative  radius. 

(7)  The  advanced  node   generation  scheme  presents  the   additional 
parameters  KRELAD,  LGTYPE,   and  extensions  to  MODEG.   Other  entries   are 
the  same  as  the  basic  generation  scheme. 

The  parameter  KRELAD  allows   one  to  specify  the  x-  and/or  y-coordinates 
of  the  current  node  by  referring  to  coordinates  of  nodes  previously 
defined.    Thus,  when  KRELAD  /  0,   the  variables   XCOORD  and/or  YCOORD 
are  assigned  the  reference  nodal  numbers  instead  of  coordinate  values. 
KRELAD  is   only  used  on  cards  with  MODEG  ■  0  or  2. 

(8)  The  parameter  LGTYPE  is  used   for  two  distinct  node   generation 
commands.   For  the  case  LGTYPE  =   1   or  2 ,   the  command  is   for  one -coordinate 
Laplace  generation,  whereas   for  the   case  LGTYPE  ■  A,   the  command  is 

for  generation  of  one-quarter  of  an  ellipse.    In  the  first  case,   one- 
coordinate  Laplace   generation  is  only  valid  for  straight  line   generation 
with  MODEG  =  2.    If  LGTYPE  -   1,   the  y-coordinate  is   specified  by  the 
straight  line  generation,  but  the  x-coordinate  is   set  by  the  Laplace 
averaging  scheme,    and  vice  versa,  if  LGTYPE  =  2. 
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The  ellipse  command  can  be  used  with  MODRG  »  2  or  3;  it  generates 
one-fourth  of  an  ellipse  between  the  node  NNP*  defined  on  the  previous 
card  and  the  node  NNP  on  the  current  card.  The  elliptical  quadrant 
is  generated  counterclockwise  with  convexity  on  the  right  when  traveling 
from  NNP*  to  NNP. 

(9)  In  addition  to  MODEG  -  0  and  2  defined  in  the  *  'basic"  scheme, 
MODEG  =  1  allows  one  to  recall  the  coordinates  of  any  previously  defined 
node,  NNP.  The  motivation  for  this  is  to  provide  a  starting  point  for 
a  new  generation  sequence  to  be  completed  on  the  following  card.  If 
the  generation  sequence  is  to  end  on  a  node  that  has  been  previously 
defined,  MODEG  ■  3  must  be  used  instead  of  MODEG  -  2.  If  MODEG  -  1 
or  3  is  specified,  XCOORD  and  YCOORD  need  not  be  input  since  these 
values  are  already  known. 

Lastly,  the  option  MODEG  =  5  permits  one  to  assign  nonsequential 
node  numbers  to  a  generated  sequence.  In  other  words,  MODEG  ■»  5  takes 
the  place  of  sequencing  the  nodal  numbers  with  the  parameter  NPINC. 
To  use  MODEG  ■  5,  first  specify  the  starting  node  as  always  done,  i.e., 
MODEG  =  0  or  1.  Then  insert  MODEG  =  5  cards  back-to-back,  entering 
each  intervening  node  number  in  the  entry  NNP.  For  the  last  node  in 
the  sequence,  use  MODEG  »  2  or  3  and  define  the  generation  desired. 

Note,  node  numbers  can  be  input  in  any  sequence;  all  nodes  un- 
specified on  nodal  cards  but  appearing  on  element  cards  will  have 
coordinates  automatically  generated  by  an  averaging  technique  called 
Laplace  generation. 
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(10)  Four  element  types  are  available  with  the  CANDE  program: 
continuum  quadrilateral,  continuum  triangle,  beam/column  (i.e.,  pipe 
element),  and  interface  element.  The  element  type  is  inherently  specified 
by  the  manner  in  which  the  element  connectivity  array,  IX(1),  IX(2) , 
IX(3),  and  IX(4) ,  is  defined: 

(a)  Quadrilaterals  -  set  IX(1),  IX(2) ,  IX(3) ,  and  IX(4)  to  the 
node  numbers  connected  to  the  element  in  counterclockwise 
order 

(b)  Triangles  -  set  IX(1)  ,  IX(2) ,  and  IX(3)  to  the  node  numbers 
connected  to  the  element  in  counterclockwise  order.  Leave  IX(4) 
=  0 

(c)  Beam/columns  -  set  IX(1)  and  IX(2)  to  node  numbers  at  the 
beam  ends,  such  that,  when  traveling  from  IX(1)  to  IX(2) ,  the 
pipe  interior  is  on  the  right.  Leave  IX(3)  «  IX(4)  ■  0 

(d)  Interface  -  set  IX(1)  and  IX(2)  to  separate  node  numbers  with 
the  same  coordinates  on  either  side  of  a  common  interface.  Set 
IX(3)  to  a  node  number  unused  by  any  other  element  (this  node  is 
associated  with  interface  forces).  Set  IX(4)  -  0  and  IX(7)  -  1 
to  distinguish  from  triangle 

These  elements  are  illustrated  in  Figure  11-10. 
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Quadrilateral 


Triangle 


pipe  interior 


Body  I 


Beam/Column 


interface 


Figure  11-10.  Types  of  CANDE  elements. 
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(11)  The  element  material  number,  IX(5) ,  identifies  groups  of 
elements  with  common  material  properties.  Material  numbers  associated 
with  one  element  type  are  unrelated  to  another  element  type,  except 
continuum  elements  (quadrilateral  and  triangle)  share  common  material 
numbers.  The  material  numbering  for  beam/column  elements  (i.e.,  pipe 
elements)  is  a  special  case  and  must  be  defined  as  follows.  Each  pipe 
element  is  given  a  separate  material  number  beginning  with  1  and  proceeding 
clockwise  around  the  pipe  in  a  connected  sequence.  This  special  numbering 
scheme  is  required  for  evaluating  the  pipe  responses. 

(12)  The  element  generation  options  require  serial  input  of  element 
numbers.  Element  numbers  that  are  skipped  between  successive  input 
cards  will  be  automatically  generated,  such  that  the  node  connectivity 
of  each  generated  element  is  determined  by  adding  INTRAL  to  the  nodes 

I,  J,  K,  and  L  of  the  previous  element  in  the  sequence.  The  material 
number  and  construction  increment  remain  as  specified  by  the  first 
element  in  the  sequence.  If  it  is  desired  to  change  the  material  number 
and/or  construction  increment  during  the  sequence  generation,  insert 
intervening  card(s)  specifying  only  NE,  IX(5) ,  and  IX(6)  of  the  element 
where  the  change  occurs,  i.e.,  do  not  input  connectivity. 

(13)  To  extend  the  generation  to  more  than  one  row  of  elements, 
simply  identify  the  number  of  rows  to  be  generated  (NUMLAY) ,  and  specify 
the  node  number  increment  from  one  row  to  the  next  INTERL  (generally 
this  would  be  the  number  of  elements  in  a  row  plus  one) . 


84 


(14)  Unless  specified  differently  on  this  card,  every  node  is 
assumed  to  be  unrestrained  with  zero  external  forces  acting  on  it  (i.e., 
IIFLG(1)  »  IIFLG(2)  =  0,  and  BIVD(1)  »  BIVD(2)  =  0.0).  To  apply  dis- 
placement constraints,  concentrated  strip  loads,  and  pressures,  repeat 
this  card  as  required.  Note,  every  node  has  two  displacement  degrees 

of  freedom.  The  boundary  condition  for  each  degree  of  freedom  is  either 
an  imposed  external  force  (IIFLG  =  0) ,  or  an  Imposed  displacement  (IIFLG 
«=  1).  The  additional  rotational  degree  of  freedom  associated  with  nodes 
of  the  beam/column  element  is  unrestrained  only  if  both  displacements 
of  the  node  are  unrestrained.  Otherwise,  a  nonrotational  end  restraint 
is  automatically  imposed. 

(15)  All  boundary  conditions  are  referenced  to  the  global  coordinate 
system.  To  reference  the  boundary  conditions  to  a  skewed  system,  define 
the  angle,  THETA. 

(16)  To  generate  a  sequence  of  identical  boundary  conditions, 
the  node  NP  is  first  in  the  sequence  and  NNP  is  the  last.  Boundary 
conditions  are  automatically  generated  for  intervening  node  numbers; 
NP  +  INCR,  NP  +  2«INCR,  etc.,  where  INCR  is  positive  for  increasing 
numbers  and  negative  for  decreasing  numbers. 

(17)  Pressure  loads  can  be  specified  over  any  connected  sequence 
of  nodes  by  specifying  a  constant  or  linearly  varying  pressure  from 
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node  NP  to  node  NNP.  Pressure  is  normal  to  each  line  segment,  and 
positive  values  point  to  left  in  traveling  from  NP  to  NNP.  For  pressure 
generation,  set  IIFLG(1)  »  IIFLG(2)  -  0. 

The  same  boundary  node  NP  may  be  specified  again  on  subsequent 
cards  under  the  following  rules:  (1)  Force  boundary  conditions  add 
into  the  BIVD  vector;  (2)  it  is  permissible  to  override  any  predefined 
force  boundary  condition  with  a  displacement  boundary  condition,  but 
not  vice  versa  (it's  good  practice  to  first  define  force  boundary  con- 
ditions, then  displacement  boundary  conditions);  (3)  boundary  conditions 
are  not  activated  until  the  load  step  equals  IA.  During  this  step  the 
value  BIVD  is  applied  to  the  system  and,  thereafter,  set  to  zero  for 
subsequent  steps. 

Every  time  a  boundary  node  is  referenced  or  generated  it  counts 
as  one  more  boundary  condition  for  computing  the  total  number  of  boundary 
conditions. 

Boundary  conditions  for  interface  element  nodes  I  and  J  represent 
displacement  degrees  of  freedom  and  are  treated  as  discussed  above. 
However,  node  K  represents  interface  forces  where  the  boundary  code 
is  defined  as: 

IIFLG(1)  ■  0,  normal  force  unknown  (typical) 
IIFLG(1)  -  1,  normal  force  specified  by  BIVD(1) 
IIFLG(2)  °  0,  shear  force  unknown  (typical) 
IIFLG(2)  -  1,  shear  force  specified  by  BIVD(2) 
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Normally  node  K  of  interface  elements  is  not  specified  in  boundary 
condition  input.  However,  if  corresponding  displacement  degrees  of 
freedom  of  nodes  I  and  J  are  both  restrained  to  zero,  then  the  corresponding 
interface  force  must  also  be  restrained  to  zero. 
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PART  III 
EXAMPLE   PROBLEMS 

The  example  problems  presented  herein  are   strictly   to  demonstrate 
data  input  and  data  output.   Accordingly,   the  problems   are   contrived 
to  illustrate  various  modeling  options.   They  are  not  necessarily  intended 
as   design  recommendations. 

The  six  examples   considered  are  identified  in  Table  III-1   by  execution 
mode,  pipe  type,   solution  level,  installation  type,    and  special  comments. 

In  subsequent  pages,  each  problem  is  presented  in  the   following 
format:    (1)    a  listing  of  all  input  data  cards  defining  the  problem, 
and   (2)    selected  CANDE  output  identifying  the  input  parameters  and 
illustrating  structural  responses   and  pipe  evaluation.   For  brevity, 
soil  response  output  is  omitted.     Unless  otherwise  noted  the  units 
of  response   data  are  in  the   inch-pound  system.    For  reference,    Table 
III-2   lists   output  symbols  with  their   definitions   and  units;   Tables 
III-3,   III-A,   III-5,   and  III-6  identify  the  design  and  evaluation  factors 
for  corrugated  steel,   corrugated  aluminum,   reinforced  concrete,   and  plastic 
pipe,   respectively. 
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Table  III-2.   CANDE  Output  Symbols 


Symbol 

Units 

Description 

ANG 

deg 

Rotation  to  coordinates  for  principal  stresses 

EPSILON-X 

in.  /in. 

Strain  in  x-coordinate  direction  (soil) 

EPSILON-Y 

in.  /in. 

Strain  in  y-coordinate  direction  (soil) 

GAMA-XY 

in.  /in. 

Engineering  shear  strain  in  x-y  coordinate 
plane  (soil) 

MOMENT 

in. -lb /in. 

Bending  moment  in  pipe  wall  section  per  unit 
length 

M- STRESS 

psi 

Outer  fiber  stress  due  to  the  bending  moment 

N-PRES 

psi 

Normal  pressure  acting  on  pipe  periphery 
(i.e. ,  soil  pressure) 

SHEAR 

lb /in. 

Resultant  shear  force  acting  on  pipe  section 

S-STRESS 

psi 

Shear  stress  on  pipe  wall  cross  section 

S-PRES 

psi 

Traction  pressure  acting  on  pipe  periphery 
(soil  shear  on  pipe) 

SIGMA- X 

psi 

Stress  in  x-coordinate  direction  (soil) 

SIGMA- Y 

psi 

Stress  in  y-coordinate  direction  (soil) 

SIGMA- 1 

psi 

Maximum  principal  stress  (soil) 

SIGMA- 2 

psi 

Minimum  principal  stress  (soil) 

TAU-XY 

psi 

Shear  stress  in  x-y  coordinate  plane  (soil) 

THRUST 

lb /in. 

Resultant  circumferential  force  in  pipe  wall 

T-STRESS 

psi 

Thrust  divided  pipe  sectional  area 

X- COORD 

in. 

Node  location  of  x-coordinate  in  undeformed 

state 

X-DISP 

in. 

Displacement  from  undeformed  position  in 
x-direction 

Y- COORD 

in. 

Node  location  of  y-coordinate  in  undeformed 
state 

Y-DISP 

in. 

Displacement  from  undeformed  position  in 
y-direction 
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Table  III-3.   Evaluation  Output  for  Corrugated  Steel 


Design  and 
Evaluation 
Items 

Definition  of  Item 

Thrust  yielding 
safety  factor 

Displacement 
collapse 
safety  factor 

Elastic  buckling 
safety  factor 

Bending  stress 
performance 

Handling 

performance 

factor 

Design  output 

Yield  stress  of  steel  divided  by  maximum 
thrust  stress  in  pipe  wall  section 

Twenty  percent  of  average  diameter  divided 
by  maximum  relative  displacement  of  pipe 

Critical  elastic  buckling  pressure  divided 
by  average  normal  pressure  experienced  by 
the  pipe 

Yield  stress  of  steel  divided  by  bending 
stress  (a  value  of  1.0  indicates  yielding 
and  is  permissible) 

Should  have  FF  >  D2/EI,  where  FF  =  0.043, 
or  0.02  for  structural  plate  corrugation 
[D  =  pipe  diameter  and  EI  =  pipe  bending 
stiffness  (elastic)] 

Wall  area  (in.2/in.) 
Moment  of  inertia  (in.^/in.) 
Section  modulus  (in.3/in.) 
List  of  suitable  corrugations 
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Table  III-4.   Evaluation  Output  for  Corrugated  Aluminum 


Design  and 
Evaluation 
Items 

Definition  of  Item 

Thrust  yielding 
safety  factor 

Yield  stress  of  aluminum  divided  by  maximum 
thrust  stress  in  pipe  wall  section 

Displacement 
collapse 
safety  factor 

Twenty  percent  of  average  diameter  divided 
by  maximum  relative  displacement  of  pipe 

Elastic  buckling 
safety  factor 

Critical  elastic  buckling  pressure  divided 
by  average  normal  pressure  experienced  by 
the  pipe 

Outer  fiber 
rupture  safety 
factor 

Strain  causing  aluminum  rupture  divided  by 
maximum  outer  fiber  strain  due  to  thrust 
and  bending 

Bending  stress 

performance 

factor 

Yield  stress  of  aluminum  divided  by  bending 
stress  (a  value  of  1.0  indicates  yielding  and 
is  permissible) 

Handling 
performance 

Should  have  FF  >  D2/EI,  where  FF  =  0.09,  or 
■  0.042  for  structural  plate  [D  =  pipe  diameter 
and  EI  =  pipe  bending  stiffness  (elastic)] 

Design  output 

Wall  area  (in.2/in.) 
Moment  of  inertia  (in.^/in.) 
Section  modulus  (in.  /in.) 
List  of  suitable  corrugations 
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Table  III-5.   Evaluation  Output  for  Reinforced  Concrete 


Design  and 
Evaluation 
Items 

Definition  of  Item 

Reinforcement 
yielding  safety 
factor 

Yield  stress  of  steel  reinforcement  divided 
by  absolute  maximum  stress  in  the  steel 
cage(s) 

Concrete  crushing 
safety  factor 

Compressive  strength  of  concrete,  f'r>    divided 
by  maximum  compressive  stress  in  outer  fibers 

Shear  cracking 
safety  factor 

Wall  shear  strength  (i.e.,  concrete  tensile 
strength)  divided  by  maximum  shear  stress  at 
axis  of  bending 

Allowable  crack- 
width  performance 
factor 

Crack  width  of  0.01  inch  divided  by  maximum 
predicted  crack  width  occurring  at  location 
of  maximum  steel  stress 

Allowable 
displacement 
performance 
factor 

o 
Allowable  pipe  deflection  =  D  /1200h,  where 

D  =  diameter  and  h  =  wall  thickness.   Factor 

=  allowable/actual 

Bow  stringing 

performance 

factor 

Tensile  strength  of  concrete  (5\/fc),  divided 
by  a  high  estimate  of  tensile  stress  between 
steel  and  concrete 

Design  output 

Required  wall  thickness  (in.) 

Required  steel  areas  for  circular  or  elliptical 
cages  (in.2/in.) 
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Table  III-6.   Evaluation  Output  for  Plastic  Pipe 


Design  and 
Evaluation 
Items 

Definition  of  Item 

Displacement 
collapse  safety 
factor 

Outer  fiber 
ultimate  stress 
safety  factor 

Elastic  buckling 
safety  factor 

Handling 

performance 

factor 

Design  output 

Twenty  percent  of  diameter  divided  by  maximum 
relative  displacement  of  pipe 

Ultimate  stress  of  plastic  material  divided 
by  predicted  maximum  outer  fiber  stress 

Critical  elastic  buckling  pressure  divided 
by  average  normal  pressure  experienced  by 
the  pipe 

Must  have  FF  >  D2/EI,  where  FF  =0.333  [D  = 
pipe  diameter  and  EI  =  pipe  handling  stiff- 
ness (elastic)] 

Required  wall  thickness  (in.) 
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Problem  1  —  Input 


INPuf  Cards  FOR  PROSLE*  NO*dER  I 

CAfiD  TtpE  J... ♦.-..  1'n  ..-..♦...  2o.-  .-.-.•*  ^<^30  <• ■  n**t. ■<*$  <<<<*,<<  50  " «,*<*«6o »*»#«#« #70#«  .,»««  ,,  8( 

CAPD  lA  DESIGN  J  STEEL      f Hi    COVER  30  f£Et,  60-I*CH  PlP£<  T-tlfT, 

CADD  IB  \                60.0 

CAPO  2B  2.0 

CAu„  1C  12C.0 

CAJO  ?C  3O.C     1000.C       O.35 
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Problem  1  -  Output 


***  PROBLEM  NUMBER   J  *** 

FILL  COVER  JO  FEET,  fcO-INCH  PIPE.  i-LIFT, 

EXECUTION  KODE  DESI 

SOLUTION  LEVEL  .....BURNS 
CULVERT  TYPE  ......  STEEL 

PIPE  PROPERTIES  ARE  AS  FOLLOWS  ... 

AVERAGE  PIPE  DIAMETER  (IN.)  .,,.,..  ,60000E»02 

YOUNGS  MODULUS  OF  PIPE  (PSI)  ,50000E*08 

POISSONS  RATIO  OP  PIPE  (-) ,3000CE*00 

YIELD  STRESS  OP  PIPE  (PSI)  ,53000E*05 

DENSITY  OF  PIPE  .(PCI)  , -0, 

MATFRIAL  CHARACTER,  NONLIN  ,...,,..     I 

NONLINsO,  IMPLIES  LINEAR 

NONLIN»l,  IMPLIES  rlELD-HlsGE 

NONLIN*2,  IMPLIES  BILINEAR  CURVE 

DESIRED  SAFETY  FACTORS  FOR  PIPE  DESIGN  ,., 

THRUST  STRESS  YIELDING  2.00600 

DEFLECTION  COLLAPSE, (0,2»D)  ..,.,,.    a. 00000 
EL*STIC  BUCKLING 2,00000 
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Problem  1  -  Output  (cont'd) 


SYSTEM  PROPERTIES  ARE  AS  FOLLOWS.,. 

DENSITY  OF  FILL  MATfRIAL  (PCF)  ....  120.000' 

NO.  OF  DATA  POINTS  ON  HALF  PIPE  ,,,  II 

NUMBER  OF  LOAD  INCREMENTS I 

SLIP/  NO  SLIP  INTERFACE  OPTION  ,,,,  NO  SLIP- 


MATERIAL  PROPERTIES  OF  SOIL  IN  VICINITY  OF 
CULVERT  AS  A  FUNCTION  OF  SOIL  HEIGHT, 


LOAD 

SOIL 

CONFINED 

L&TERIAL 

YOUNGS 

POISSONS 

STEP 

HEIGHT 

MODULUS 

STRESS 

MODULUS 

RATIO 

(FEET) 

(PSI) 

RATIO 

(PSI) 

30.00         1600.94  ,50      1000,00        ,350 
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Problem  1  -  Output  (cont'd) 


DESIGN  PESULTS  AT  2  ITERATIONS 

**«  REQUIRED  THRUST  ARE*  .........     .05691 

*•♦  REQUIRED  MOM,  OR  INERTIA  , 00277 

**♦  REQUIRED  SECTION  MODULUS  ,,,.,     .00*00 


THE  FOLLOWING  STEEL  CORRUGATED  SECTIONS  MEET  THE 
A80VE  REQUIREMENTS  WITH  MIM^M  area 


CORRUGATION 

GAGE     THRUST 

ARE* 

MOM, 

INERTIA 

3/2  x 

I/O 

(NO  ADEQUATE 

SECTION) 

2  x 

1/2 

12 

,11900 

,0035a 

8/3  X 

J/2 

12 

,11300 

,00302 

1  X 

1 

18 

,05925 

,00689 

6  X 

2 

12 

. 129  6,6 

,06001 

THE  CORRUGATION  with  minimum  AREA  AMD 
MOMENT  OF  INERTIA  IS    3X1    GAGE.  18 
AN  ANALYSIS  OF  THIS  CORRUGATION  FOLLO-S. 
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Problem  1  -  Output  (cont'd) 


STRUCTURAL  RESPONSES  OF  CULVERT 


10*0  STEP  ■  1 


X-COORD, 
Y-C00R0, 

X-0ISP, 
Y-OISR, 

N-PRES, 
S-PRES, 

MOMENT 
M. STRESS 

THRUST 
T-STRESS 

SHEAR 
S-STPESS 

•  ,00 
30,00 

.105E-16 
•.04UE+00 

•.1O6E+02 
.03UE-13 

,363E*03 
,276E*05 

,53PE*03 
.91 0E*0fl 

-.636E-13 
-.  i  o  1  E-l 1 

«,27 
26.53 

.115E-01 
•,fl26E*00 

•.206E+02 
•,7UflE*0l 

,293E»03          • 
.223EOS 

,578E*03 
,976E»0a 

,1C3E*02 
,2fl2£*03 

17.63 
2«,27 

.925E-01 
•,266E*00 

-,23!E*02 
-,120E*02 

,UCE*03 

,637E*0O           « 

,679E»03 
,115E*05 

,232E»02 
,3«1E*0S 

2«,27 
17.63 

,2U7E«00 
•,106E«00 

-,263E»02 
-,120E*02 

«,U6E*03           « 
•,661E*0O 

,80aE«03 
•,136E*05 

.232E+02 
,391E*0S 

28,53 

«.27 

,ao5E*oo 

•.186E-01 

-,289E*02 
•,7fl<JE*01 

-,29OE»03           « 
•,227E*05 

>,905E*0S 
■,15JE*05 

,lO3E*02 

,2fl2E*03 

30.00 
•  ,00 

,O71E*00 
.762E-16 

«.298E*02 
.135E-12 

•,368E«03 

-,260E*05           < 

•,9aaE*03 
>.159E*05 

-.260E-12 
•.038E-U 

26,53 

-9,27 

,a05E*00 
.186E-01 

-,26OE»02 

,7auE+01 

•,29OE*03 
-,227E»05 

».905E*03 
•,153E*05 

-,103E»02 
-,2A2E*03 

2a, 27 
•17,63 

,2U7E«00 
,106E*00 

-,263E*02 

,120E*02 

•,U6E*03 
•,881E*0O 

►,8Q8E*03 
•,13bE*05 

-,232E*02 
-,391E*03 

17.63 
-2a, 27 

.925E-01 
,266E*00 

-,231E*02 
,120E*02 

,M0E*03          < 
.837E40U 

•,679E»03 
»,115E*05 

•,23?f»02 

•,391E*03 

9.27 
•26.53 

.115E-01 
.826E+00 

•,206E*02 
,7U«E*01 

,293E«03 
,223E*05 

»,578E*03 
.,976E*04 

-,1«3E»02 
-,2«2E»03 

,00 
•30,00 

•.057E-16 

,U9<!E*00 

•  ,1<S6E*02 
.137E-12 

,363E*03 
,276E*05 

■,539E*0S 

»,<>10E*0tt 

-.263E-12 
•.aa3E-ll 
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Problem  1  —  Output  (cont'd) 


CALCULATED  SAFETY  FACTORS  AT  STEP   i 

THRUST  STRESS  ( YIELD-STRESS/*AX. STRESS)  ■  2,0? 

DISPLACEMENT  ...,..,,.  CDIA^ETER^O, 2/hax-DISP.)  ■  12, U 

ELASTIC  BUCKLING  (P-CRITICAUP-AVERAGE)  «  6,90 

PERFORMANCE  FACTORS 

BENDING  STRESS  ,,..,,  ( YIELO-STRESS/MAX-STRESS)  ■  1,18 

HANDLING  REQUIREMENT  (FF.EI/D»«23  ■  2,09 
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Problem  2  —  Input 


input  cards  for  problem  nuvqer    ? 

CARD  TtPE  l...*...U.  ....... 20.... «... 30....*... 4  0 50. ...... .60. ...... .70 ....... .80 

C4RD  ]A  ANALYS  2    STEEL     OVERBURDEN  LO&OINC-.  3  X  \    CORRUGATION  )2-GAGE. 

CAPD  18  2  60.0 

C*«0  2ti  0-13    0»0l545    o*''?798 

Ca°D  1C  E^Ba     EMr^^em  mesi-  *1th  SOFT  eACKPACMNG. 

CAHD  ?C  1        -I  32.5  25.0 

CA"D  3C         )o 

Ca"d  Ifj  1     1  iNSjTu  aiNr^R) 

Ca*-D  ?0  1000.0       0.35 

C*RO  lD  ?.  i  BEDDING  «  iNSlTf 

CA*D  ?0  1C0C.O       0.35 

CAPO  10  3     1  FILL  »  INSITL 

CARD  20  1000.0       0.35 

CA^D    lD  L         4  1  PACKING    «    INSITL/2 

C*HD    ?D  500.0  0.35 


101 


Problem  2  —  Output 


***  PROBLEM  NUMBER   2  «*• 

OVERBURDEN  LOADING,  3  X  I  CORRUGATION  12-GAGE, 

EXECUTION  MODE  .....ANAL 
SOLUTION  LEVEL  .....f.E.AUTO 
CULVERT  TYPE  ......  STEEL 

PIPE  PROPERTIES  ARE  AS  FOLLO*S  ... 

AVERAGE  PIPE  OIAMETER  (IN.)  ,,  ,60000002 

YOUNGS  MODULUS  OF  PI?E  (PSI)  ......  ,30000C«0» 

POISSONS  RATIO  OP  PIPE  (•)  ,  ,30000E*00 

YIELD  STRESS  OP  PIPE  (PSI)  .33000EOS 

DENSITY  OP  PIPE  (PCI)  -0, 

MATERIAL  CHARACTER,  NONLIN  ,,     2 

NOnLINjO,  IMPLIES  LINEAR 
N0NLIN=1,  IMPLIES  YIELD-HINGE 
NONLINaJ,  IMPLIES  BILINEAR  CURVE. 

MODULUS  OF  UPPER  BI-CURVE  .........    »0, 

SECTION  PROPERTIES  OF  PIPE  ,., 

THRUST  AREA  (IN««2/IN) ,13000 

MOM,  OF  INERTIA  (IN««0/IN)  ,,,     ,01505 

SECTION  MODULUS  (IN**3/IN)  ,,,, 0279$ 
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Problem  2  -  Output  (cont'd) 


•  *  BEGIN  GENERATION  OF  CANNED  MESH  •  # 

THE  DATA  TO  BE  RUN  IS  ENTITLED 

EMBANKMENT  *ESH  WITH  SOFT  BACKPACKING. 

TYPE  OF  MESH— —.———. ...........  EMBANKMENT 

PLOTTING  DATA  SAVED • — — —  -0 

PRINT  SOIL  RESPONSES t 

PRINT  CONTROL  FOR  PREP  OUTPUT----------  -0 

NUMBER  OF  CONSTRUCTION  INCREMENTS--  —  -  1 

PIPE  DIAMETE«  RATIO 1,00 

LIVE  LOAD  PRESSURE  ON  LAST  INCR, 25.00 

NUMBER  OF  BACKPACKING  ELEMENTS-  — •  —  --  10 

WIDTH  OF  BACKPACKING  LAYER — --       5.00 

SOIL  HEIGHT  ABOVE  SPRINGLIKE  (FEET)--—  32,50 

MESH  HEIGHT  ABOVE  SPRINGLIKE  (FEET)--—  10,00 

SOIL  DENSITY  ABOVE  MESH  (PCF) 0,00 

PRESSURE  OF  TRUNCATED  SOIL  (PSD— — —      0,00 

IDENTIFICATION  OF  MATERIAL  ZONE  WITH  MATERIAL  NUMBERS, 

HATERIAL-ZONE    MATERIAL  NO, 

INSITU  1 

BEDDING  2 

BACKFILL  3 

BACKPACK  4 
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Problem  2  -  Output  (cont'd) 


*  *  BEGIN  PREP  OF  FINITE  ELEMENT  INPUT  •  » 

THE  DATA  TO  BE  RUN  15  ENTITLED 

EMBANKMENT  KESH  WITH  SOFT  BACKPACKING, 

NUMBER  CF  CONSTRUCTION  INCREMENTS------  1 

PRINT  CONTROL  FOR  PREP  OUTPUT--  — -----  -0 

INPUT  DATA  CHECK--  —  -----  —  ---  —  —  ----  0 

PLOT  TAPE  GENERATION---  —  —  —  —  —  —  —  -0 

ENTIRE  FINITE  ELEMENT  RESULTS  OUTPUT-—  1 

THE  NUMBER  OF  NODES  IS- - — -•  110 

THE  NUMBER  OF  ELEMENTS  IS--------------  96 

THE  NUMBER  OF  BOUNDARY  CONDITIONS  IS---  OS 
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Problem  2  -  Output  (cont'd) 


******  FOR  MESH  GENERATION  OP  *****  * 


EMBANKMENT  KESM  WITH  SOFT  BACKPACKING. 


THE  NUMBER  OK  DATA  ERRORS  IS-----  —  —  -  0 

THE  NUMBER  OF  SOIL  MATERIALS  IS  — 4 

THE  NUMBER  OF  PIPE  MATERIALS  IS--------  10 

THE  NUMBER  OF  INTERFACE  MATERIALS  IS—  0 

THE  BANO  WIDTH  IS Zb 

*******  MESH  DATA  HAS  BEEN  SAVED  *  *  * 


AN  OROERED  LIST  OF  BEAM  ELFMENTS  AROUND  PI°E. 

ORDERING  IS  CLOCKWISE  FPOK  TOP  TO  BOTTOM  OR  LEFT  TO  RITE, 

12305678*        10 
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Problem  2  -  Output  (cont'd) 


MATERIAL  CKARAC.R1ZATI0N  FOR  SOILS. 


PROPERTIES  FOR  MATERIAL   i    ****»«♦*  INSITU  (LINEAR) 

DENSITY  i    -0, 

YOUNGS  MODULUS*  .1000EOU 

POISSONS  RATIO*  ,3500E*00 

CONFINED  MOD,«  .  1605E*00 

LATERAL  COEFF.a  ,5385E*00 

PROPERTIES  FOR  MATERIAL   2   *«»***««BEDDING  ■  INSITU 

DENSITY  a    -0, 

YOU'JGS  MODULUS*  ,1000E*04 

POISSONS  RATIOS  ,3500E*00 

CONFINED  MOD,*  ,UCSE*0U 

LATERAL  COEFF.a  ,5385E*00 

PROPERTIES  FOR  MATERIAL   3    *»*»*•♦»    FILL  ■  INSITU 

DENSITY  i    »0, 

YOUNGS  MODULUS*  .1000S*0« 

POISSONS  PATIO*  ,3500E*00 

CONFINED  MOD,*  ,1605E*0a 

LATERAL  COEFF.a  .5385E+00 

PROPERTIES  FOR  MATERIAL   «    *«**»«»*PACKING  »  INSITU/2 

DENSITY  *    -0. 

YOUNGS  MODULUS*  ,5000E*03 

POISSONS  RATIO*  .3500E+00 

CONFINED  MC0.»  ,8025E*03 

LATERAL  COEFF.a  ,538SE*00 
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Problem  2  -  Output  (cont'd) 


STRUCTURAL  RESPONSES  OP  CULVERT 


LOAD  STEP  ■  1 


X-COORO, 
Y-COORD, 

y-oisp, 

Y-OISP, 

N»PRES. 
S-PRES, 

HO^ENT 
M-STRESS 

THRUST 
T-STRESS 

S«E*R 
S-STRESS 

0.00 
30,00 

o. 

-,232E*0l 

-,219E*02 
0. 

,809E*03 
.289E+05 

•,502E*03 
•,017E*00 

0. 
0. 

9.27 
28,53 

,90"E-02 
•.225E+01 

-,222E*02 
-.515E-.01 

,650E*03 
,232E*0S 

•,5*1E*03 

•,uo7E*0a 

,301E*02 
.232E+03 

17.63 
2*1.27 

.8S3E-01 
-,210E*01 

«.225E*02 
-,816E*01 

,2O3E*03 
,869E»00 

-,657E»03 
>.50SE»0a 

,O7OE»02 
,3b8E*03 

20.27 
17.63 

,233E*00 
-,195E*01 

-,2flOE*02 
•,8|flE*01 

-,2O9E*03 
•,889E»0O 

>.7O9E»03 
•,S7bE*0O 

,O77E»02 
.367E+03 

28.53 

9.27 

,38UE*00 
-,187E*0l 

-.2O6E+02 
-,099E*01 

-,652E*03 
-,233E*05 

».82flE»03 
■,630E*0« 

,296E»02 
,227E*03 

30,00 

0.00 

,OU7E*00 
-,186E*01 

-,209E*02 
,101E»00 

-,80«E*03 
-.287E+05 

».851E*03 
•,655E*0a 

-,3a3E*00 
-,2baE»01 

28,53 
•  9,28 

,362E*00 
•,185E*01 

-,206E»02 
,500E*01 

-,6U5E*0.3 
-,231E*05 

».82?£»03 
•,632E*0a 

-,29RE*02 
-,230E»03 

20,27 
•17,60 

,231E*00 
-,i77E*01 

-,239E*02 

,8JOE*01 

-,2fl2E»03 

-,86lE»0U 

•,7fl7E»03 
■,575F»0a 

-,O7PE*02 
-,36<>E«03 

17,63 

-20,27 

.8OJE-01 
•,162E*01 

-,230E«02 

,8ioe»oi 

,25UE*03 
,909E»0a 

»,655E*03 

•,50aE*Oa 

-,872E»02 
-.563E03 

*.27 
•26,53 

.883E-02 
•.1O7E+01 

-,221E»02 
,O92E*01 

,6uoE»03    « 
,230E*05 

>,561E»03 
»,uu7E*0O 

-,266f*02 
-,220E»03 

0.00 
-30,00 

0. 

•,1ooe*oi 

-,217E*02 
0, 

,79JE»03 
,283E»05 

»,5U0E*03 
■,0J8E*0o 

0. 
0, 

107 


Problem  2  -  Output  (cont'd) 


NONLINEAR  CHARACTERISTICS,    LOAD  STEP  «   1 


OE 

INNFR-PIBER 

OUTEP»EISER 

STRAIN    PATIO 

FRACTION    OF 

STRAIN 

STRAIN 

MAX" 

•TOYIE10 

HALL 

YIELDED 

1 

,750"?E-05 

-.10035E-02 

l,00?ao 

,00102 

2 

.56937E-03 

»,8ao<J9E»03 

,83065 

0,00000 

3 

,! 10U9E-03 

-,<U687E«03 

,M6U6 

0.00000 

a 

•  ,<jauu  1E>03 

,9a?25E-0O 

,ua3<»6 

0,00000 

5 

•.89B98E-03 

.51U57E-03 

,89508 

0,00000 

6 

•, 1 0700E-0? 

.6727UE-03 

1.06898 

,039fc2 

7 

•.891U8E-03 

.5079JE-03 

,89059 

0,00000 

6 

-.U3618E-03 

.87593E-CU 

,C3570 

0,00000 

0 

.1230SE-03 

-.U2863E-03 

.U2820 

0,00000 

10 

,56278E»03 

-.83387E-03 

,83303 

0,00000 

It 

.73067E-03 

«,98UuuE»03 

.<>83a5 

0.00000 
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Problem  2  —  Output  (cont'd) 


CALCULATED  SAFETY  FACTORS  AT  STEP   ! 

THRUST  STRESS  ,, ( Y IELD-STRESS/mAX-STRESS)  «  5,00 

DISPLACEMENT  ,,,,  (DIAMETER*0,2/MAX-OISP,)  *  13.15 

ELASTIC  BUCKLING  ..,,...  (P-CRITICAL/P-AVERAGE)  e  10,35 

PERFORMANCE  FACTORS 

BENDING  STRESS  ,  ( YIELD-STRESS/MAX-STRESS )  a  1.10 

HANDLING  REQUIREMENT  ,..,.,,, (FF*EI/D»*2)  «  5,57 


109 


Problem  3-  Input 


IKPyT  CAROS  FOP-  PROBLEM  N(j"HER   3 

CARD  TyPe       1... ♦• ..IT ........ 20. ....... 30.  • ...... 40. t».».»»50  ••••♦•••*<>• •♦.♦••♦70. .♦.♦..• 86 

FILL  COVER  30  PF.ET  IN  3  LIFTS,  FlXfO  WftLL  OESIGN, 
6«o  3 

<-0000» 
C.01        FIX 


CAPO  1A 

OESI6N  1  cON'CR£ 

CARD  jB 

&C«0 

CARU  2b 

CAPD  3B 

CARO  ,C 

120T   .  3 

CARD  ?C 

10.0 

CARO  ?C 

?°-° 

CARD  ^C 

3fl.o    1 

7CO.0 

0.35 

900.0 

0.35 

loco.o 

0.35 
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Problem  3  -  Output 


*♦*  PROBLEM  NUMBER   3  •#* 

FILL  COVER  SO  FEET  IN  3  LIFTS,  FIXED  W»LL  DESIGN, 

EXECUTION  MODE  ,..,.DESI 

SOLUTION  LEVEL  BURNS 

CULVERT  TYPE  ,..,..  CONCRETE 
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Problem  3  -  Output  (cont'd) 


PIPE  PROPERTIES  ARE  AS  FOLLOWS  ,,, 
(UNITS  ARE  INCH-POUND  SYSTEM  ) 

NOMINAL  PIPE  DIAMETER  .,  60,0000 

CONCRETE  COMPRESSIVE  STRENGTH  ,,,,,  0000,0000 

CONCRETE  ELASTIC  MODULUS  3830253,5127 

CONCRETE  POISSON  RATIO  .,.,.,,,,,..  .1700 

STEEL  YIELD  STRENGTH  , f,,,  60000,0000 

STEEL  ELASTIC  MODULUS  29000000, 0000 

STEEL  POISSON  RATIO  ,,,., ,,  ,3000 

NONLINEAR  CODE  (1,2, OR  3)  ,  3 

CONC,  CRACKING  STRAIN  (1,2,3)  0,000000 

CONC,  YIELDING  STRAIN  (2,3)  ,000507 

CONC,  CRUSHING  STRAIN  (2,3)  ,002000 

STEEL  YIELDING  STRAIN  (3)  ,001883 

DESIRED.  SAFETY  FACTORS  FOR  P^E  DESIGN  ,,, 

STFEL  YIELDING  1,6000 

CONCRETE  CRUSHING 2.0000 

CONCRETE  SHEAR  FAILURE  ,.,,..  2,0000 

ALLOWABLE  CRACK  WIDTH  (IN)  .,,,,,,  .0100 

8HAPE  OF  REBAR  CAGE(S)  CIPC 

MIN,  COVER  TO  CENTER  OF  REBAR  1,2500 

DESIGN  GOAL  FOR  WALL  THICKNESS  ,,,,  FIX 

FIXED  WALL  THICKNESS  6,0000 

OUTER-TO-INNER  STEEL  PATIO  ..,,.,..  ,7500 
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Problem  3  —  Output  (cont'd) 


SYSTEM  PROPERTIES  ARE  AS  FOLLOWS,,. 

DENSITY  OF  FILL  MATERIAL  (PCF)  ,,.,  120,000 

NO.  OF  DATA  POINTS  ON  HALF  PIPE  ...  11 

NUMBER  OF  LOAD  INCREMENTS 3 

SLIP.  NO  SLIP  INTERFACE  OPTION  ....  NO  SLIP 


MATERIAL  PROPERTIES  OF  SOIL  IN  VICINJTY  OF 
CULVERT  AS  A  FUNCTION  OF  SOIL  HEIGHT. 

LOAD         SOIL  CONFINED  LATERAL  YOUNGS    POISSONS 

STEP       HEIGHT  MODULUS  STRESS  MODULUS      RATIO 

(FEET)  (PSI)  RATIO  (PSI) 

1  10.00  1123.06  .50  700,00        .350 

2  20,00  latu.oa  ,50  900.00       ,350 

3  30,00  160U,9a  ,50  1000,00       ,350 
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Problem  3  —  Output  (cont'd) 


DESIGN  RESULTS  »,   5  ITERATIONS 

#*»  HALL  THICKNESS  , 6.00006 

•  •*  INNER  CAGE  STEEL  AREA  COR  ELLI)     .03<J7a 

•  »*  OUTER  CAGE  STEEL  AREA 02<»8J 
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Problem  3  -  Output  (cont'd) 


STRUCTURAL  RESPONSE  OF  CULVERT   FOR  LOAD  INCREMENT 


COORDINATES,  DISPLACEMENTS  AND  CRACK  DEPTHS  ARE  IN  INCHES 
PRESSURES  ARE  IN  LB/IN**2 
MOMENTS  ARE  IN  IN,*LB/IN, 
THRUST  AND  SHEAR  ARE  IN  LB/IN, 


NPPT 


10 


11 


X-COORD, 
Y-COORD, 

X-DISP, 
Y-OISP. 

N-PPES, 
S-PRES, 

HOHENT 
THRUST 

SWEAR 
CRACK  DEPTH 

-,00 
35.00 

.57095E-17 
-.R5766E-01 

-,5i595E*02 
,5a272E-15 

,01902E*Oa 
-,52976E*03 

-.87965E-12 
,a2958E*01 

10,20 
31,38 

.20618E-02 
-.82676E-01 

-,50323E»02 
-,56727E*0t 

,53821E*0fl 
-,585T5E*03 

,15073E»03 
,a2022E»01 

IR.ao 
26,70 

.17517E-01 
-.51616E-01 

-,2b993E*02 
•,95O22E*0l 

,12bb5E»0a 
-,73227E»03 

,2aj89E»03 
,58639E*0l 

26,70 
19, ao 

,a72aoE-01 
•.20697E-01 

-,2287feE*02 
-,95022E»01 

-,13e85E*08 
-,9j3aoE«03 

,2a589E*05 
,3b058E*01 

31,58 
10,20 

.77552E-01 
-.3733SE-02 

-,19586E*02 
-,58727E*01 

-,38fcaiE*0u 
-,10599E»0a 

,1S073E*03 
,ai800E»01 

35.00 

-,00 

.90357E-01 
.17500E-16 

-,J827aE*02- 
.106U9E-12 

-,a2722E»0a 
-.11 15«E*0a 

-.27332E-I1 
,a2519E*01 

31,30 

•10,20 

.77532E-01 
.37333E-02 

•,195e6E*02 
,58727E40i 

-,3«6aiE*oa 
-,10599E*Ca 

-.15073E+03 
,ai600E*01 

26,70 
•19. ao 

,O72a0E«01 
.20697E-01 

-,2287bE»02 
,95022E*01 

-,13a85E»08 
-,R13aoE*03 

-,2aJ89E*03 

,3b05BE*01 

19, ao 
•26,70 

.17517E-01 
.51616E-01 

-,2b993E*02 
,95022E*01 

,12bb5E*0a 
-,73227E»03 

-,'2a389E»03 
,58b39E*01 

10.20 
-31.36 

.20618E-02 
.82676E-01 

-,30323E»02 
,S8727E«0l 

,33621E*0O 

-,58573E*03 

-,!5073E*03 
,a2022E*0l 

,00 
•35,00 

-.11658E-16 
.95766E-01 

-,31595E*02 
.10771E-12 

,ei902E*oa 
-,52976E»03 

-,27babE-ll 
,82O5aE*01 
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Problem  3  -  Output  (cont'd) 


STRESSES  IN  CULVtRT  HALL   (PSI)   FOR  LOAD  INCREMENT 


ELLIP.  OR 

PT 

INNER  CAGE 

OUTER  CAGE 

CONCRETE 

SKEAP 

STEEL 

STEEL 

COMPRESSION 

STRESS 

1 

,2277aE*05 

•,33997E*0a 

•>.1579UE*0« 

-.21228E-12 

2 

,1779GE*05 

-,33013E*0<J 

•,13"?3E*0a 

,J6979E»0? 

3 

,«7t J9E+0U 

«,27383E*0U 

•.66R02E+03 

,709<J<iE*02 

a 

•  t32207E<-0U 

,66309E*OU 

-,83u9feE*03 

,2«»«11E*02 

5 

•,4SP8PE*0U 

,2255<5E  +  05 

-,17362E*0U 

,57772E*01 

6 

•■,U7U53E*00 

,28599E*05 

-,20215E»0U 

-,C62bfcE-13 

7 

•t43epsE*oa 

.2255OE+05 

-,17362E*0U 

-,57772E*01 

6 

-,32207E*0« 

,66309E*0U 

»,83t96E*03 

•  .29i-UE*02 

9 

,<I71  3<?r  *o« 

•t?7383E*0a 

-,fc6<'f|2E*03 

•,7099aE*P2 

10 

,17790E*0S 

-.33013EOU 

•,13«23E*0O 

•,36979E»02 

11 

,2277«E«05 

-,33997EfOO 

-.1579uE«0fl 

•.66716E-12 
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Problem  3  -  Output  (cont'd) 


CALCULATED  SAFETY  FACTORS  FOR  LOAO  INCREMENT      3 

STEEL  YIELD  STRESS  /  MAX.  STEEL  STRESS  2,0*8 

CONCRETE  STRENGTH  /  max.  CO"PPESSIVE  STRESS  ....  1.97"» 

WALL  SMEAR  CAPACITY  /  MAX,  SMEAR  ..,.,..,.,,,,.,  0.050 
PERFORMANCE  FACTORS 

0.01  INCH  /  MAX,  CRACK  WIDTH  .,,,,, ...•,..  .'TO 

ALLOCABLE  DISPLACEMENT(LUM)  /  max.  DISP,  ,,,.,..  2. fell 

TENSILE  STRENGTH  /  STRESS  FROM  90*  STRING  , ,«18 
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Problem  4—  Input 


INPUT  cAR0S  fOK  Pf-CBLfM  NtM°ER   4 

CAPO  TtPE       l«'«*,"lO«..«*«..20««««*«««3o«».«*«.«4g....«.»«5o»»..4«t.6o.««»*.«.7o«»««*«».8o 

CARD  1A  ANALY5  2  PLaSTI     fill  cOvfR  »  20  FEET,  48-INCH  P^E.  FPICTION  FACE. 

CAoO  IB  48.0 

CAon  26  0.50 

CA«D  1C  homo    '  INCLUDE  INTEPFacE  ELEMENTS  FCR  FRICTION  *  0.1  SLIP 

SSS  '  22-°      120-c 

CA»D  ID 
CAPD  20 
CAPO  ,0 
CAOD  20 
CARD  iD 
CARD  2° 
CA°D  10 

card  ?n 
CapD  id 
CARD  20 
CAPO  ID 
CARD  2D 
CARD  10 
CARD  20 
CAPD  ID 
CA"D  20 
CARD  jD 
CArO  ?D 
CAPD  10 
CARD  2D 
CAPO  ID 
CAPD  20 
CARD  ID 
CAPD  2D  90<n        OM 


1 

-1 

22.0 

p 

1        1 

120. C 

HOWCGENEOUS    SOIL 

IGOO'O 

0'35 

1            6 

-90.  ft 

0.1 

2           6 

"72-0 

OM 

3           6 

-5a. n 

O'l 

4            6 

-36.-5 

0.1 

5           6 

-18. C 

0.1 

6           6 

0.00 

°M 

7           6 

le.o 

0.1 

8           6 

36.0 

a.  I 

9            6 

54. 0 

0.1 

10           6 

72.0 

0.1 

11           6 
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Problem  4  —  Output 


*«•  PROBLEM  NUMBER   0  *** 

rill  COVER  »  20  FEET,  08-INCH  PIPE,  FRICTION  FACE, 

EXECUTION  KOOE  ANAL 

SOLUTION  LEVEL  ,.,..F,E, AUTO 
CULVERT  TYPE  ......  PLASTIC 

PIPE  PROPERTIES  ARE  AS  FOLLOWS  ... 

AVERAGE  PIPE  0IA"ETER  (IN.)  OB, 00 

YOUNGS  MODULUS  OF  PIPE  (PSI)  1600000.00 

POISSONS  RATIO  OF  PIPE  (-)  ,30 

YIELD  STRESS  OF  PIPE  (PS!)  25000.00 

DENSITY  OF  PIPE  (PCI) »0.00 

SECTION  PROPERTIES  OF  PIPE  ... 

THICKNESS  OF  PIPEWALL  (IN) ,50000 

THRUST  AREA  (IN»«2/IN)  ., ,50000 

MOM,  OF  INERTIA  (IN»»a/IN)  ,.,,..,,  ,01002 

SECTION  MODULUS  (IN«»3/IN)  ,, ,  ,00167 
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Problem  4  —  Output  (cont'd) 


*  *  BEGIN  GENERATION  OF  CANNED  MESH  *  * 

THE  DATA  TO  BE  RUN  IS  ENTITLED 

INCLUDE  INTERFACE  ELEMENTS  FOR  FRICTION  i  0,1 

TYPE  OF  MESH—. HOMOGENOUS 

PLOTTING  DATA  SAVED"-  — .««  -0 

PRINT  SOIL  RESPONSES-.- j 

PRINT  CONTROL  FOR  PRE°  OUTPUT.-...— «--  -0 

NUMBER  OF  CONSTRUCTION  INCREMENTS-  —  ...  1 
PIPE  DIAMETER  RATIO  —  -..-».--  — .......      1.00 

LIVE  LOAD  PRESSURE  ON  L'ST  INCR. .-  —  -0.00 

SOIL  HEIGHT  ABOVE  SPRINGLINE  (FEET)--—  22,00 
HESH  HEIGHT  ABOVE  SPRINCLINE  (FEET)-—      8,06 

SOIL  DENSITY  ABOVE  MESH  (PCF) 120,00 

PRESSURE  OF  TRUNCATED  SOIL  (PSI)— — ••  11,60 

IDENTIFICATION  OF  KATERIAL  ZONE  WITH  MATERIAL  NUMBERS, 


MATERIAL-ZONE  MATERIAL  NO, 

INSITU  I 

BEDDING  1 

BACKFILL  1 
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Problem  4  —  Output  (cont'd) 


*  *  BEGIN  PREP  OF  FINITE  ELEMENT  INPUT  *  « 

THE  OAT*  TO  BE  RUN  IS  ENTITLEO 

INCLUDE  INTERFACE  ELEMENTS  FOR  FRICTION  u    0,1 

NUMBER  OF  CONSTRUCTION  INCREMENTS------  1 

PRINT  CONTROL  FOR  PREP  OUTPUT---------  -0 

INPUT  DATA  CHECK---------------  -  ---  0 

PLOT  TAPE  GENERATION 0 

ENTIRE  FINITE  ELEMENT  RESULTS  OUTPUT---  1 

THE  NUM8ER  OF  NODES  IS-----------------  132 

THE  NUMBER  OF  ELEMENTS  IS —  — .-  107 

THE  NUMBER  OF  BOUNDARY  CONDITIONS  IS---  ««? 
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Problem  4  -  Output  (cont'd) 


******  FOR  MESH  GENERATION  OP  *  *  *  *  •  * 

INCLUDE  INTERFACE  ELEMENTS  FOR  FRICTION  «  0.1 

THE  NUMBER  OF  DATA  ERRORS  IS--  —  —  -—    0 

THE  NUMBER  OF  SOIL  MATERIALS  IS •     1 

THE  NUMBER  OF  PIPE  MATERIALS  IS  —  —  -    10 
THE  NUMBER  OF  INTERFACE  MATERIALS  IS  —  -    Ji 
THE  BAND  WIDTH  IS  —  ---  —  -------  —  —    tb 

*******  MESH  OATA  HAS  BEEN  SAVED  *••••* 


AN  ORDERED  LIST  OF  BEAM  ELEMENTS  AROUND  PIPE, 

ORDERING  IS  CLOCKWISE  FROM  TOR  TO  BOTTOM  OR  LEFT  TO  RITE, 

123956769    10 
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Problem  4  —  Output  (cont'd) 


MATERIAL  CHARACTERIZATION  FOR  SOILS, 

PROPERTIES  FOR  MATERIAL   1    «»•••»•»     HOMOGENEOUS  SOIL 

DENSITY  ■      ,12000E*03 

YOUNGS  MOOULUSc  ,1000E*0« 

POISSONS  RATIO"  ,3500E*00 

CONFINED  ►'CD, a  ,160S£*0O 

LATERAL  COEFF.s  ,5365E*00 

INTERFACE  ELEMENT  HATERI AL-GROUP  PROPERTIES 
MAT,  NO,    NORMAL-ANGLE   COEF-FPICTION   TENSILE-RUPTURE 


t 

-90,00000 

,10000 

,00000 

2 

•72,00000 

,10000 

.00000 

3 

•SU. 00000 

,10000 

,00000 

a 

-36,00000 

,10000 

,00000 

5 

•16,00000 

,10000 

.00000 

6 

0,00000 

.10000 

,00000 

7 

18,00000 

,10000 

,00000 

8 

36,00000 

,10000 

,00000 

9 

5fl.00C00 

,10000 

,00000 

10 

72.00000 

,10000 

.00000 

11 

90,00000 

,10000 

,00000 
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Problem  4  —  Output  (cont'd) 


STRUCTURAL  RESPONSES  OF  CULVERT 


LOAD  STEP  «  1 


X-COORD, 
V-COORD, 

X-DISP. 

Y-DISP, 

N-PRES, 
S-PRES, 

HOMEKT 
h-STRESS 

THRUST 
T-STRESS 

SHEAR 
S-STRESS 

0,00 
2U.25 

0. 
".156E+C1 

0 

,166E*02 

,301E*02 
,723E*03 

>.39JE»03 
-,783E*03 

0. 
0. 

7.U9 

23,06 

■350E-02 
-.152E+01 

m 
• 

,lb8E«02 
,168E*01 

,237E*02 
,569E*03 

>,A03E*03 
•,806E»03 

,10CE*01 
.280EO1 

1«,25 
19,62 

.521E-01 
-,IU1E*01 

- 

,172E*02 
,172E*01 

,885E*01 
,212E*03 

..fll7E*03 

',83aE*03 

,200E*01 
,O01E»01 

19,62 
1«,25 

,150E*00 

-.131E+01 

• 

,176E*02 
,178E*01 

-,670E»01 
«,161E*03 

•,O3lE*03 
-,8b2E*03 

,222E*01 

,fl83E»0l 

23,06 
7,«9 

,251E*00 
-,125E*01 

* 

,162E*02 
,182E*01 

-,2a?£*02 
-,595E*03 

»,ec5E*03 
>,69iE*03 

,157€*01 
,313E»01 

21,25 
0,00 

.300E+00 
-,12UE*01 

- 

,  1  6RE-»02 
,168E*01 

-,305E*02 
-,731E*03 

■,a60E*03 
•,920E*03 

,fl53E-01 
.9C7E-01 

23,06 

-7,50 

.259E+00 
-.123E+01 

- 

188E+02 
.188E+01 

-,25SE*02 
-,611E*03 

■,OfcOE*03 
',919E*03 

-,131E*01 

-,262E«01 

19.62 
•1«,26 

,156E*00 
-.117E+01 

- 

,162E*02 
,182E*01 

-,10t>E*02 
-,253E*03 

■,885E*03 
».890E*03 

-,23af*01 
'.its  SE  ♦  0  1 

11.25 
•19,62 

.505E-01 
•.107E+01 

• 

,178E*02 
,178E*01 

,100f*02 
,2UOE*03 

>.830£O3 
»,860E*03 

•,23aE*01 
-,«67E»01 

7.U9 
-23.06 

.363E-02 
-,959E*00 

- 

,173E*02 
,173E*01 

,2fl9E*02 
,597E*03 

»,ai6E*03 
'.832E*03 

-,139E»01 
-,278E*01 

0.00 
•20.25 

0. 
-,913E*00 

0 

,171E*02 

,311E*02 
,706E*03 

-,008E»03 
'.808E*03 

0. 

0. 

124 


Problem  4  —  Output  (cont'd) 


CALCULATED  SAFETY  FACTORS  AT  STEP   1 

DISPLACEMENT  ., (DIAMETER*0.2/HAX-DISP.)  »  15,00 

OUTER-FIBER  STRESS  ...  (ULTIMATE  /  MAXIMUM)      «  15.01 

ELASTIC  BUCKLING  (P-CRITICA|./P-AVER AGE)  ■  3. TO 

PERFORMANCE  FACTORS 

HANOELING (FF»E»I/D**2)  «  2,36 


125 


Problem  5  —  Input 


INPUT  CiROS  FOR  PROBLEM  NtMgER   5 

C  A  f<  D  TtPE  j«..*»«.)r;'»««*«««2o*,******30*,,,*,*,*0********5o,,,**,,"'>0******,,^0,,*,**,,3( 

CAPO  1A  DESI    2  ALUmIN  1«ENCh  OEPTh  9  FEET.  TRENCH  WIDTH  7  FEET,  OVERFILL  15  FEET 

CAPO  ]8  2    £,0.0 

CARD  2B  2.5        4.0        2.0        2.0 

CADD    1C  TPEN      *00IFY  L£V£L  2  FOp  CCNCirNTfJ4TEO  STRIP  L0AD  »  400lBs/IN  (TOTaL>      MOd 

card  ?c  ii       -l             i5.n    i2o.o 

CArO  3C  9.0        7.0 

CARD  4C  0     0     1 

CARD  7C  I03                        -200."               1 

CApll  ]C  x             j             INSITU  <LINFAR» 

CARD  20  1U00.0    0.3s 

CARD  ID  2     1            BEODING  »  InSItu 

CARD  ?D  lCOO.O    0.35 

CARD    ID  L         3           1       120.0           TRENCH   FILL    =    1*00 

CARD    20  1500.0         0.3S 
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Problem  5  —  Output 


•*•  PR0BLEM  NUMBER   5  *** 

TRENCH  DEPTH  9  FEET,  TRENCH  WIDTH  7  FEET,  OVERFILL  15  FEET 

EXECUTION  MODE  ,,...DESI 

SOLUTION  LEVEL  F.E.AUTO 

CULVERT  TYPE  ALUMINUM 

PIPE  PROPERTIES  ARE  AS  FOLLOWS  ,., 

AVERAGE  PIPE  DIAMETER  (IN.) ,  ,60000E»02 

YOUNGS  MODULUS  OF  PIPE  (PS!)  ,,,,.,  ,10200E*CS 

POISSONS  RATIO  OF  PIPE  CO  ..,,.,..  ,330001*00 

YIELD  STRESS  OF  PIPE  (PSl')  ..,.,,,,  ,2flOO0E«05 

RUPTURE  STRAIN  OF  PIPE  (IN/IN)  ,,,,  .50000E-01 

DENSITY  O'  PIPE  (PCI) -0, 

MATERIAL  CHARACTER,  N0"LIN  2 

NONLINtO,  IMPLIES  LINEAR 
NONLINsl,  IMPLIES  YIELO-MINGE 
NONLIN:2,  IMPLIES  BILINEAR  CURVE 

MODULUS  OF  UPPER  BI-CURVE ,60000E»06 

DESIRED  SAFETY  FACTORS  FOR  PIPE  DESIGN  ... 

THRUST  8TRESS  YIELDING  .,.,.,.,,.,,  2,50000 

DEFLECTION  COLL*pSE, (0,2«D)  a, 00000 

ELASTIC  BUCKLING  2,00000 

OUTER  FIBER  RUPTURE  ,,,,  2,00000 
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Problem  5  —  Output  (cont'd) 


*  *  BEGIN  GENERATION  OF  CANNED  MES*  *  * 

THE  DATA  TO  BE  RUN  IS  EN'ITLEO 

MODIFY  LEVEL  2  FOR  CONCENTRATED  STRIP  LOAD  a  O00L5S/IN  (TOTAL) 

TYPE  OF  KESH- ....................... ..  TRENCH 

PLOTTING  DATA  SAVED°»-----------------  J 

PRINT  SOIL  RESPONSES---- 1 

PRINT  CONTROL  FOR  PREP  OUTPUT-———  -0 

NUMBER  OF  CONSTRUCTION  INCREMENTS------  1 

PIPE  DIAMETER  RATIO - I. 00 

LIVE  LOAD  PRESSURE  ON  LAST  INCR. ------  -0,00 

DEPTH  OF  TRENCH———————  9,00 

WIDTH  OF  TRENCH........................  7,00 

SOIL  HEIGHT  ABOVE  SPRINGLINE  (FEET)-  —  -  21,50 

HESH  HEIGHT  ABOVE  SPRINGLINE  (FEET)----  6.50 

SOIL  OENSITY  ABOVE  MESH  (PCF) 120,00 

PRESSURE  OF  TRUNCATEO  SOIL  (PSI) 12,50 

IDENTIFICATION  OF  MATERIAL  ZONE  *ITH  MATERIAL  NUMBERS, 


HATERIAL-ZONE     MATERIAL  NO, 


INSITU 

BEDDING 
TRENCH 
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Problem  5  —  Output  (cont'd) 


•  *  BEGIN  PREP  OF  FINITE  ELEMENT  INPUT  «  ♦ 

THE  DATA  TO  BE  RUN  15  ENTITLED 

MODIFY  LEVEL  I    FOR  CONCENTRATED  STRIP  LOAD  «  OOOLBS/IN  (TOTAL) 

NUMBER  OF  CONSTRUCTION  INCREMENTS------  1 

PRINT  CONTROL  FOR  PREP  OUTPUT------ ---.  -0 

INPUT  DATA  CHECK--------- ...  —  .»  0 

PLOT  TAPE  GENERATION-.----.-..,.-.-.--.  J 

CNTIPE  FINITE  ELE"ENT  RESULTS  OUTPUT---  1 

THE  NUMBER  OF  NODES  IS 110 

THE  NUM8ER  OF  ELEMENTS  IS— -•----•--•--  96 

THE  NUMBER  OF  BOUNDARY  CONDITIONS  IS---  C5 
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Problem  5  —  Output  (cont'd) 


•  •  *  CHANGES  TD  STANDAPD  LEVEL  2  MESH  »  *  *  ♦ 

a  NUMBER  OF  NODf"^,  TO  BE  CHANGED  .--  —  .-       0* 

*  NUMBER  OF  ELECTS  TO  BE  CHANGED  — -       0* 

*  ADDITIONAL  BOUNDARY  CONDITIONS  ••-  —  •       1*, 

a******************   •*•* 


♦••ADDITIONAL  BOUNDARY  CONDITIONS,, .FORCES  ■  LBS   ,  DISPLACEMENTS  »  INCHES... 

BOUNDARY  LOAD       X-FORCE      OR       T-FORCE      OR       X«Y  ROTATION 

NODE  STEP      X-DISPLACEMENT      Y-DISPLACEHENT  DEGREES 

103         I     r  «   »0.  f   «   ».2000E*05       -0, 
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Problem  5  —  Output  (cont'd) 


FOR  MESH  GENERATION  OF  *  *  •  *  *  * 


MODIFY  LEVEL  2  FOR  CONCENTRATED  STRIP  LOAD  c  flOOLBS/IN  (TOTAL) 

THE  NUMBER  OF  DATA  ERRORS  IS  —  —  —  —  —     0 

THE  NUMBER  OF  SOIL  MATE»IALS  IS --•     3 

THE  NUMBER  OF  PIPE  MATERIALS  IS———    10 
THE  NUMBER  OF  INTERFACE  MATERIALS  IS—     0 
THE  BAND  WIDTH  IS—————————    26 

*******  KESH  DATA  HAS  BEEN  SAVED  •*•*••• 


AN  ORDERED  LIST  OF  PEA*  ELEMENTS  AROUND  PIPE, 

ORDERING  IS  CLOCKWISE  FPOM  TOP  TO  BOTTOM  OR  LEFT  TO  RITE, 

123056769    10 
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Problem  5  -  Output  (cont'd) 


MATERIAL  CHiRACltRIZATION  FOR  SOILS, 

PROPERTIES  FOR  MATERIAL   1    *»**»**»  INSITU  CLINEAR) 

DENSITY  «    -0, 

YOUNGS  wODULUS=  ,1000E*00 

POISSONS  »ATIOs  .3500E+00 

CONFINED  *OD,»  ,1605E*0a 

LATERAL  COEFF.b  ,5365E*00 

PROPERTIES  FOR  MATERIAL   2    ***#»*««  BEDDING  »  INSITU 

DENSITY  ■    -0, 

YOUNGS  MODULUS*  ,1000E*Oa 

POISSONS  PATIO*  ,3500E*00 

CONFINED  MOD,:  ,lfc05E*0U 

LATERAL  COEFF.s  .538SE+00 

PROPERTIES  FOR  MATERIAL   5    ***••*«*  TRENCH  FILL  *  1500 

DENSITY  «      ,12000E*OJ 

YOUNGS  MODULUS!  ,150C£*0« 

POISSONS  RATIO*  ,350CE*00 

CONFINED  MOD, a  ,2<i07E*0<i 

LATERAL  COEFF.a  ,53B5E»00 
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Problem  5  —  Output  (cont'd) 


DESIGN  RESULTS  AT  2  ITERATIONS 

«**  REQUIRED  THRUST  AREA  ,,.,,,,,.  ,06965 

***  REOUIREO  MOM,  OR  INERTIA  ,,,,,  .00392 

*«*  REQUIRED  SECTION  MODULUS  .....  ,01067 


THE  FOLLOWING  ALUMINUM  CC»RUG4TED  SECTIONS  MEET  THE 
ABOVE  REQUIREMENTS  *<ITH  HJNJHum  *R£a 

CORRUGATION     GAGE     THRUST  AREA  MO*.  INERTIA 

6/3  X   1/2       10         ,1«533         ,00fl53 
3X1  16  ,07016  .00865 

6X1  14  ,08080  ,01060 

6X2  6  ,20«08  ,09616 

9X5/2         0  .11700  ,09100 


THE  CORRUGATION  KITH  minimum  AREA  AND 
MOMENT  OF  INEPT1A  IS    3X1    GAGE,  16 
AN  ANALYSIS  OF  THIS  CORRUGATION  FOLLO«S, 
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Problem  5  —  Output  (cont'd) 


STRUCTURAL  RESPONSES  OF  CULVERT 


LOAD  STEP  *  1 


X-COOPD, 
Y-COORD, 

X-DISP. 
Y-DISP, 

N-PPES, 
S-PRES, 

ypME\T 
K-STRESS 

THRUST 
T-STRESS 

SMEAR 
S-STRESS 

0,00 
30,00 

0. 

-,176E*01 

-.109E+02 
0. 

,182E*03 
,lilE*05 

■,2«5E*03 
'.385E*0a 

0. 
0. 

9.27 
2B.53 

.817E-02 
-.169E+01 

-.118E+02 
-,96?E*01 

,125E*03 
,767E*0O 

•,337E*03 

•ta5aE*oa 

,867E*01 
,117E»03 

17,63 
20.27 

,8a5E-01 
-.153E401 

-,1U0E*02 
-,13<>E*02 

,157E*02 
,U5E*0a 

-,«52E*03 
•«610E»0a 

,10CE*02 
.135EOJ 

2«, 27 
17,63 

,220E*00 
-.138E+01 

-,1R0E*02 
-.117E+02 

-,631E*C2 
-,3P6E*0« 

».577E*03 

»,778E*0a 

,7fe«E»01 

,10aE«03 

2*. 53 
9.27 

.3fc«E*00 
-,131E*01 

-,213E*02 
-,557E*01 

-,126E*03 
-,76<JE*0« 

•,661E»03 
»,891E*0a 

,337E*01 
,O55E«02 

30,00 
0,00 

.363E+03 
-.130E+01 

-,217E»02 
,30<?E*01 

-,126E*03 
-,77uE*0O 

■,671E*03 
■,<?C5E*0a 

-,301E*0! 
-,a05E*02 

28,53 

•<»,2e 

.322E+00 
-.128E+01 

-,210E*02 
,O5aE*01 

-,6"»2F*02 
-,a2«;E*0a 

».608F*03 
•,820E*08 

-,?38E*01 
-,32tE«02 

2«,27 
-17,6a 

.205E+00 
-,121E*01 

-,155E*02 
,108E*02 

-,817E*02 

-,500E*0a 

■,510E*0J 

-,688E»0U 

-,560E*01 
-,761E»02 

17.63 
-2«,27 

.750E-01 
«.108E*01 

-.1O0E+02 
,101E*02 

,3fc7E*02 
,225E*0O 

■,aoeE»03 
•,S51E*0a 

-.1OE402 
-,1«7E*03 

9.27 
-26,53 

.620E-02 
-.SJjE+OO 

-,117E*02 
,573E*01 

,123E»03 

,751E*0O 

-,3J1E»03 

■,aa&E*0a 

-,571E»01 
-,76«E*02 

0,00 
-30,00 

0. 
-.672E+00 

-,106E*02 
0. 

,iaaE*03 
,660E*00 

■.2<?8E*03 
•,602E*oa 

0. 
0. 
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Problem  5  —  Output  (cont'd) 


NONLINEAR  CHARACTERISTICS,    LOAD  STEP 

DE 

INNER-FIBER 

OUTER-FIBER 

STRAIN 

STRAIN 

1 

,63<J7SE»03 

-.13066E-02 

2 

.27285E-03 

-,10b6<?E-02 

3 

-.432S5E-03 

-,632b9E-03 

0 

-.10168E-02 

-.3U237E-03 

5 

•  .U505E-02 

-,10bb3c»03 

6 

•tia660E-02 

-,ll<J«9E-03 

7 

-.10863E-02 

-,3«il8E-0S 

6 

•.10376E-02 

-,lb"15E-03 

9 

•.28UbbE-0S 

-.67735E-03 

10 

,26bb3E-03 

-.10U55E-02 

11 

,U1737E'03 

-.1120UE-02 

STRAIN  BATIO 
MAX-TO-YIELD 


FRACTION  OF 
WALL  YIELDED 


,62319 

0.00000 

,50887 

0,00000 

.30175 

O.COOOO 

, UgUQO 

0,00000 

,b9[78 

o.oocoo 

,69<519 

0,00000 

.51811 

0,00000 

.U9U59 

0,00000 

,32306 

0,00000 

,U<J8b6 

0,00000 

.53*36 

0.00000 
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Problem  5  —  Output  (cont'd) 


CALCULATED  SAFETY  FACTORS  AT  STEP   1 

THRUST  STRESS  ( Y1ELD-STRESS/HAX-STRESS)  «  2,65 

DISPLACEMENT  , (DIAMETER*0.2/HAX-D!SP,)  *  13. 5S 

ELASTIC  BUCKLING  (P-CRI T ICAL/P-A VERAGE )  «  7,57 

STRAIN  RUPTURE  ....  (RUPTURE-STR AIN/*AX«STPA IN)  «  30,11 

PERFORMANCE  FACTORS 

BENDING  STRESS  ,,,,.,  (YlELD-STRESS/MAX-STRESS)  ■  2.16 

HANDLING  REQUIREMENT {FF«EI/D««2)  «  2,21 
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Problem  6-  Input 


CARD  TtPE 

CARD  lA 
CAPO  >B 
CARD  2B 
CAPO  3B 
CA°D  1C 
C*pD  2C 
CApO  3C 
CAPO  3C 
CAoO  3C 
CAHO  *C 
CARD  AC 
CAPO  AC 
CAwO  4C 
CAPO  4C 
CAPO  AC 
CApO  «C 
CARD  AC 
CARD  «C 
CARD  4C 
CAB!)  5C 
CAPO  5C 
CARD  5C 
CAH0  5c 
CARD  5C 
CA^D  5C 
CAPU  ^c 
C*pO  ]C 
CAPU  7U 


INPUT  CARUS  f°W  PP'JBLEM  NUKRe"   6 

,.♦... 20 3O 40 50....*...  60....*...  70....* 


ANAL.YS  3  CONCPE 


6C.0 
Anon. 
0,03666     C.C275 
PREP*     D  LOAD  -fCSf 


0  LOAD  TEST,  fo-INCH  OIA.,  CLASS  3t  *AlL  A. 
5.0  ClP-C    3 


6 
1 
6 
Jl 
1 
2 
3 

4 

5 
6 

7 

8 
9 


1 
2 

3 
4 

5 

fc 
7 
8 
9 
10 


.  .  .80 
1011 


0.0 

30.0 

O.C 


1.0 


11 

in 

3C 

11 

.0 

.0 

-3(1 

.0 

1 

2 

3 
4 

5 

6 
7 

8 

9 

10 

-7?« 

0 

-7«;< 

c 

•75 

p 

-75- 

n 

►75 

0 

• 

'75 

f 

1.0 

7 


-30.0 
•30.0 


0.0 


0*0 
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Problem  6  -  Output 


•**  PROBLEM  NUMBER   6  •** 

0  LOAD  TEST,  60-INCH  DIA,,  CLASS  3,  WALL  A, 

EXECUTION  MODE  anal 

SOLUTION  LEVEL  F.E.USER 

CULVERT  TYPE  CONCRETE 
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Problem  6  -  Output  (cont'd) 


PIPE  PROPERTIES  APE  AS  FOLLOWS  ... 
(UNITS  APE  INCH-POUND  SYSTEM  ) 

NOMINAL  PIPE  OIAKETEP  .............  60.0000 

CONCRETE  COMPRESSIVE  STRENGTH  ,,,,,  0000. 0000 

CONCRETE  ELASTIC  MODULUS  383a253.5l27 

CONCRETE  POISSON  PATIO ,.  .1700 

STEEL  YIELD  STRENGTH  .....,...,...,  flOOOO.OOOo 

STEEL  ELASTIC  "ODULUS  ,   2*000000.0000 

STEEL  POISSON  BATIO  , .5000 

NONLINEAR  CODE  (1,2, OR  3)  .........  3 

CONC,  CRACKING  STRAIN  (1,2,3)  0,000000 

CONC.  YIELDING  STRAIN  (2,3) ,000507 

CONC,  CRUSHING  STRAIN  (2,3)  .002000 

STEEL  YIELOING  STRAIN  (3)  ,001255 

WALL  THICKNESS  ....................  5.0000 

SHAPE  OF  REBAR  CAGES  ...,.,,.,,...,  CIRC 

INNER  CAGE  STEEL  AREA  ,0367 

OUTER  CAGE  STEEL  AREA  ,,,,, ,0275 

COVER  TO  CENTER  OF  INNER  PEBAR  ....  1,0000 

COVER  TO  CENTER  OF  OUTER  REBAR  ,,,,  1,0000 
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Problem  6  —  Output  (cont'd) 


*  *    BEGIN  PREP  OF  FINITE  ELE*ENT  INPUT  *  * 

THE  DATA  TO  BE  RUN  IS  ENTITLED 
D  LOaD  TEST 

NU*!3ER  OF  CONSTRUCTION  INCREMENTS------  6 

PRINT  CONTROL  FOR  PREP  OUTPUT——-—  fl 

INPUT  DATA  CHECK----- ------------------  -0 

PLOT  TAPE  GENERATION-  — ---- .0 

ENTIRE  FINITE  ELE"ENT  RESULTS  OUTPUT.-.  »0 

THE  NUMBER  OF  NODES  is------------- ---  11 

THE  NUMBER  OF  ELEMENTS  I S------ --------  10 

THE  NUMBER  OF  BOUNDARY  CONDITION'S  IS---  7 
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Problem  6  —  Output  (cont'd) 


******  FOR  KESH  GELATION  OF  *****  * 

D  LOAD  TEST 

THE  NU«BE5  OF  DATA  ERRORS  IS-----------     0 

THE  NUMBER  OF  SOIL  MATERIALS  IS 0 

THE  NUMBER  OF  PIPE  MATERIALS  IS---— ---    10 
THE  NUMBER  OF  INTERFACE  MATERIALS  IS---     0 

THE  BAND  WIDTH  is---------- --....-     a 

*******  KESH  DATA  HAS  BEEN  SAVED  *•••••• 


AN  ORDERED  LIST  OF  BEAM  ELEMENTS  AROUND  PIPE. 

ORDERING  IS  CLOCKWISE  FROM  TOP  TO  BOTTOM  OR  LEFT  TO  RITE, 

1  23«5t>769        10 
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Problem  6  —  Output  (cont'd) 


MATERIAL  CHARACTERIZATION  FOR  SOILS, 


PROPERTIES  FOR  MATERIAL   1    ******** 

DENSITY  ■    -0. 

YOUNGS  MODULUS*  0, 

POISSONS  RATIO:  0, 

CONFINED  MOD. a  0, 

LATERAL  COEFF.s  0, 

KARNIVC  MORE  MATERIAS  WERE  INpUT  THAN  CEFINEO     1     0 
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Problem  6  -  Output  (cont'd) 


STRUCTURAL  RESPONSE  OF  CULVERT   FOR  LOAD  INCREMENT 


COORDINATES,  DISPLACEMENTS  AND  CRACK  DEPTHS  ARE  IN  INCHES 
PRESSURES  ARE  IN  LB/IN**2 
MOMENTS  ARE  IN  IN,*LB/IN, 
THRUST  AND  SHEAR  ARE  IN  LB/IN, 


NPPT 

X-COORD. 
Y-COORD, 

X»DISP. 
Y-DISP, 

N-PRES, 
S-PRES. 

MOHENT 
THRUST 

SHEAR 
CRACK  DEPTH 

1 

0.00 
30,00 

0. 

-.56609E-01 

-,15981E*02 
0. 

.loeoRE^oo 

,320taE-09 

0. 
,37592E*01 

2 

9.27 
26.53 

.7J8U3E-03 
-.516S7E-01 

-.2369RE-U 

•,215b5E-10 

,78<?«>2E*03 
-,22891E*02 

,70fl51E*02 
,37112E»01 

3 

17.4.5 

20.27 

.56718E-02 
-,02121E»01 

•,21608E«10 
-.J6802E-10 

,1623<»E*03 

»,O35u:E«02 

,5»R2<JE*02 
,3322flE*01 

a 

20.27 
17.63 

,103UAE»01 
-.33361E-01 

,22aeuE.10 
•.06256E-J3 

-,335S5E»03 
-.59  929E+02 

.O35«lE*02 
,36fl97E>01 

5 

28.53 
0.27 

.22762E-0! 
-,28»73E-01 

-,  1 039PE-1 0 
-,26<583E»1S 

•,650Q7E*03 
-,70U51E*02 

,226<?1E*02 
,37«23EOl 

6 

30.00 
0,00 

.26215E-01 
«.2832SE«0i 

-.60396E-11 
•.71220E-11 

•,7650«5E»03 
-,7«077E*02 

-,a<>u5aE»n 
,37566E*0t 

7 

28,53 

•9,27 

.22762E-01 
-.27676E-01 

•.1O736E-10 
•.15767E-1I 

-,65097E*03 
-,70<i5lE*02 

-.22891E*02 
.37O23E+01 

0 

20.27 
•17.63 

.1O308E-01 
•.23268E-01 

•.76089E-11 
•.73727E-12 

•,33538E*03 
-,59929E»02 

•,O35aiE*02 
,36a«J7E*0l 

9 

17.63 
•20,27 

.56718E-02 
-.IO526E-01 

-.85603E-11 
,599a«E«ll 

,!6239E«03 
-fu35ulE»02 

-,5'»<»29E*02 
,3322«E*0t 

10 

'.27 
•26,53 

.73803E-03 
•.07627E-02 

-.50O97E-12 
.59715E-11 

,78962E*03 
•,22891E»02 

•,7C051E*02 
,37112E*01 

11 

0,00 
•30,00 

0. 
0. 

-,15961E*02 
0. 

.10SaOE*0a 
-.17553E-09 

0. 
,375R2E»01 
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Problem  6  —  Output  (cont'd) 


STRESSES  IN  CULVERT  HALL   (PSD   FOR  LOAD  INCREMENT 


ELLIP,    OR 

PT 

INNER    CAGE 

OUTER    CAGE 

CONCRETE 

SHEAR 

STEEL 

STEEL 

COMPRESSION 

STRESS 

1 

,11<I00E*05 

•  ,<?o«65E»03 

«.fc35!7E*03 

0. 

2 

,5<J567E*0U 

-,63UJ<JE*03 

•,35080E*C3 

,1<?852E*02 

3 

,101?<5E*00 

-,2O5(;<5E*03 

- ,  o  e  fc  a  5>  E  *  o  2 

.18C«feE02 

a 

»,U09SeE*03 

,30O7t,E  +  0u 

-.  JOSSOE  +  OS 

,l0P5«EOl 

5 

-,5P875E*03 

,626U0E*0O 

-.355<55E*03 

,3<55l  3E»00 

6 

-,6aR5JE*03 

,735StE*0tt 

- .  « 1 1 11E*03 

-.  P06P5E-13 

7 

-,58e75E*03 

,626«0E*CU 

-.  355<!5S*03 

•  ,3"513EO0 

e 

-.aosseE+os 

,3CR7feE*0u 

-,1'?55<'E*03 

-,JOS5aE*01 

9 

.101 29E*oa 

-,295act*03 

-,<50bt5E*02 

-,16C<J6E»C2 

10 

,59567E*0U 

-,fc3u3<;E  +  03 

-,35C6CE«03 

-,l«»e52E»02 

li 

,HO00E*05 

».<»<J«eSEt03 

-,63517E»03 

0. 
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Problem  6  —  Output  (cont'd) 


CALCULATED  SAFETY  FACTORS  FOR  LOAD  INCREMENT      1 

STEEL  YIELD  STRESS  /  HAX,  STEEL  STRESS  , ,,,  3.509 

CONCRETE  STRENGTH  /  MAX,  COMPRESSIVE  STRESS  ,,,,  6,298 

WALL  SHEAR  CAPACITY  /  MAX,  SHEAR  ,, ...,.*.  15.929 

PERFORMANCE  FACTORS 

0,01  INCH  /  MAX,  CRACK  WICTH  ,,,,,, ......  3.809 

ALLOCABLE  DISPLACEMENT(LUM)  /  MAX.  DIS",  .,.,,,,  10.591 

TENSILE  STRENGTH  /  STRESS  FPOM  60**  STRING  ,.,,.,  2.059 
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Problem  6  —  Output  (cont'd) 


STRUCTURAL  RESPONSE  OF  CULVERT   FOR  LOAD  INCREMENT 


COORDINATES,  DISPLACEMENTS  AND  CRACK  DEPTHS  ARE  IN  INCHES 
PRESSURES  ARE  IN  LB/IN*»2 
HOHFNTS  ARE  IN  IN,«LB/IN, 
THRUST  AND  SHEAR  ARE  IN  LB/IN, 


KPPT 


10 


11 


x-coord, 

VCOORD, 

X-DISP, 
Y-DISP, 

N-PRES. 
S-PRES, 

HOHENT                          SHEAR 
THRUST           CRACK    DEPTH 

0,00 
30,00 

0. 
•,1130CE*00 

-,310fe2E*02 
0, 

.206*oE*Ott           o 
.10O77E-08 

[375O3E*01 

0.27 
28,53 

,ta761E-02 
-,10386E»00 

-.13130E-10 
-.U2371E-10 

,15783E»0O 
-,fl5782E*02 

,lfl0R0E*03 
,37113E*0i 

17.63 
2U.27 

,113flOE-0i 
•.8O32CE-01 

•.17253E-10 
-.25710E-I0 

.32J87E03              , 

-,87082E*02              , 

11R86E*03 
.33372EO! 

20.27 
17,63 

.28722E-01 
-.66770E-01 

.27O75E-10 
•.O2603E-10 

•,67167E*03 
•.11086E403              , 

,67082r*02 
36501E»01 

26.53 
R. 27 

.05563E-01 
-.57O9U-01 

-.1O35UE-10 
-.78378E-10 

»,13109E*0a 
-,ia0O0E*03              , 

,a5782E*02 
,37a2OE*0l 

30,00 
0.00 

.52U73E-01 
•.566ORE-01 

-.135O2E-10 

-.25872E-10 

-,15311E*00              , 
•,10815E»03              , 

,1552QE-0O 
37568E*01 

28,53 
•0,27 

.U5563E-01 
-.55309E-01 

-.27185E-10 
-.15065E-10 

•,1310OE*0fi 
-,iaO<90Et03              , 

C5782E*02 
370JOE*01 

2a. 27 
•17,63 

.28722E-01 
•.O6618E-01 

•.173ORE-J0 
-,a7363£-li 

-,67167E»03 

•  ,I108iEO3             , 

e?082E*02 
36S01EO1 

17.63 

•2a, 27 

.1I3U9E-01 
-.20073E-01 

-,ia370E-10 
.171U8E-10 

,32387E»03 
-,87C82E«02              , 

I1986F*03 
33372E*01 

0.27 
•28,53 

,ia761E»02 
-.O5220E-02 

,aa3<UE-12 
.I01U6E-10 

,1578JE»0O           -, 
-,«5782E*02             , 

IflOOCE^OJ 
3T113E»01 

0,00 
•30,00 

0. 
0, 

-,S1062E*02 
0. 

,2Ofe8OE*0O   ,        0, 
-.53751E-0O              , 

375O3E*01 

146 


Problem  6  -  Output  (cont'd) 


STRESSES  IN  CUU..KT  WALL   (PSD   FOR  LOAO  INCREMENT 


ELLIP.    OP 

PT 

INNER    CAGE 

OUTER    CAGE 

CONCOETE 

SHEAR 

STEEL 

STEEL           COMPRESSION! 

STRESS 

,227RaE*05 

-,lR887E*0a 

,126R<»E*0a 

0, 

,11<>37E*C5 

«.1267RE*0u 

,7C11°E*03 

,3R70flE*02 

,2C610E*oa 

-.58G«<?E»03 

,16*.b7r*0J 

,35<JbCE*02 

-.81<)01E*03 

,6202«E*P« 

,3<>!"5E*03 

,2167flE*?l 

-,  M777E*0a 

,!25J9E*05 

,71?38E»03 

,7?«77E»00 

•.  129<»lE*0a 

. !  u  7  2  u  E  ♦  C  5 

.822*>6E*03 

,2u5fe8E-l  1 

-,11777E*0U 

,12S3<>E*05 

,7123*E*93 

-,78'977E*Ce 

-,eiR01E*CJ 

,6202aE*0a 

,3<;:<i5E*03 

•,2U7«E»0! 

,20M0c  +  0a 

•  ,58aU<JE*03 

,:eit>7E*o3 

-,35960E*92 

10 

,HR07E*05 

».12fc7<>E*0a 

,70UQE*03 

-,3«70aE*02 

11 

,227R«E*05 

•.1R887E+04 

,12bR«E*04 

0. 
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Problem  6  -  Output  (cont'd) 


CALCULATED  SAFETY  FACTORS  FOR  LOAD  INCREMENT      2 

STEEL  YIELD  STRESS  /  MAX,  STEEL  STRESS  1.755 

CONCRETE  STRENGTH  /  Mix,  COMPRESSIVE  STRESS  ,,,,  3,150 

**LL  SHEAR  CAPACITY  /  MAX,  SMEAR  ,,,.,.,,,.,,,,,  7,965 
PERFORMANCE  FACTORS 

0,01  INCH  /  MAX,  CRACK  WlOTH  ,,,,,.,,,,,.,.,,,,,  1.380 

ALLOWABLE  DISPLACEMENTCLUM)  /  MAX,  DISP,  ,  5,291 

TENSILE  STRENGTH  /  STRESS  FPOM  BOw  STRING  ,  1,230 
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Problem  6  -  Output  (cont'd) 


STRUCTURAL  RESPONSE  OF  CULVERT   FOR  LOAD  INCREMENT 


COORDINATES,  DISPLACEMENTS  AND  CRACK  DEPTHS  ARE  IN  INCHES 
PRESSURES  ARE  IN  LB/IN»»2 
MOMENTS  APE  IN  IN,»LB/IN, 
THRUST  AND  SHEAB  ARE  IN  LB/IN, 


NPPT 

X-COORD, 

X-DISP, 

N.PRES, 

HOMENT 

SHEAR 

V-COORD. 

Y-DISP, 

S-PRES, 

THRUST     C> 

*ACK  OEPTrl 

1 

0,00 

0. 

•,07903E^02 

,aa52RE*C0     0 

30,00 

•,17015E*00 

0, 

.17753E-06 

,37593E«01 

2 

".27 

,2213*E-02 

.76053E-12 

.23670E+0O 

,2U35E*03 

28.53 

-,15587E*00 

-,6eC30E-10 

-,b8673E*02 

,37113E*01 

3 

17.63 

.17027E-01 

•.26268E-10 

,O8533E*03 

,I79TCE*03 

20.27 

•,12653E*00 

-.O7591E-10 

-,13062E*03 

,33021E*01 

a 

20,27 

.O3096E-01 

.231O7E-10 

-,I0060E*0O 

,13062E*03 

17,63 

-,10C20E*00 

-.769COE-10 

-,17<>79E»03 

,3fe502E*01 

5 

28.53" 

.6836UE-01 

•.15978E-10 

-,t96b?E*00 

,6?673E»02 

".27 

-.8702UE-01 

•,630ulE-10 

-,21135E*03 

,37U25E»01 

6 

30,00 

.787J2E-01 

•.195I3E-10 

-,22971E*0O 

,2»a67E-09 

0.00 

•,6507UE»01 

-.60922E-10 

-,22223E»03 

,3756«E*01 

r 

28.53 

.6836UE-01 

-.O3685E-I0 

-,19668E*0a     • 

68673E*02 

•  9,27 

-.8312UE-01 

-.23160E-10 

-,21135E»03 

37O25E*01 

6 

20,27 

,030<56E»01 

-.16552E-10 

•,10080E»0a 

13C62E»03 

•17.63 

-,699U9E»01 

.10252E-10 

",1797«E*C3      , 

36502E*01 

9 

17,63 

.17027E-01 

".76177E-11 

,0853JE»03 

17<>79E*03 

•20,27 

-.O3618E-01 

.19JB9E-10 

•,130b2E»03      , 

33021E»01 

10 

9.27 

.22138E-02 

•.70976E-12 

,23670E»0O 

21135E»03 

•26,53 

-.1O281E-01 

.13998E-10 

-,68673E*02 

37113E»01 

11 

0,00 

0. 

-,079«3E*02 

,aa52<5E»0O    0 

•30,00 

0, 

0. 

-.7266RE-09      , 

37593E*01 
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Problem  6  -  Output  (cont'd) 


STRESSES  IN  CULVERT  WALL   (PSI)   FOR  LOAD  INCREMENT 


NPPT 


1 
2 

3 

o 
S 
6 

7 

e 

10 

>i 


ELLIP.  OR 
ISNER  CAGE 
STEEL 


.3018 

.1T85 

.3109 

-.122? 

'.  1  766 

•,19«8 

-.1766 

-.1228 

.3100 

,  1765 

,3U18 


PE  +  05 
7E*05 
8E  +  0U 

5E*oa 

feE  +  Oa 
6E*0U 
fct  +  00 
5E  +  0a 
6E»0U 
7E*05 
8E  +  05 


OUTER  CAGE 
STEEL 


-.2982 

■,1<?01 
■,8735 
,0307 
.1881 
.2200 
,1»81 
,0307 
'.8735 
'.1O01 
',2082 


5E40a 
3E»0a 
CEt-03 
5E  +  0a 
3E*05 
1E*05 
3E*05 
5E*0u 
0E*03 
JE-Oa 
5E  +  08 


CONCRETE 
CO"PRESS10N 


•.1O0U6E 
-.10516E 

-.2727JE 
•.587JJE 
-,10fc88E 
•,123«SE 
-.1C686E 
•,587J3E 
-.272735 
-,1051tE 
-,190a6E 


♦  oa 

♦  Oa 

♦  03 

♦  03 
*C« 

♦  oa 

♦  oa 

♦  03 

♦  03 

♦  oa 

♦  oa 


SHEAR 
STRESS 

0, 

,5OS5bE^02 

,5387aE*02 

,32ao2E+01 

,1!8JCE»01 

,flb52frE-ll 

-,lieasE+01 

-,32ao?E.3i 

-,5337aE»02 

-,5O556E^02 

0. 
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Problem  6  -  Output  (cont'd) 


CALCULATED  SAFETY  FACTORS  FOR  LOAD  INCREMENT     3 

STEEL  YIELD  STRESS  /  MAX,  STEEL  STRESS  J. 170 

CONCRETE  STRENGTH  /  MAX,  COMPRESSIVE  STRESS  ....  2.100 

WALL  SHEAR  CAPACITY  /  MAX.  SHEAR  ,.,, ,,.,  5.310 

PERFORMANCE  FACTORS 

0,01  INCH  /  HAX,  CRACK  WIDTH  , „„,„  ,640 

ALLOWABLE  DISPLACEMENT(LUM)  /  MAX.  DISP 3.526 

TENSILE  STRENGTH  /  STRESS  FROM  SOW  STRING  ,  ,820 
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00057171, 


FHWA 


R&D 


